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Superfluid has been realized in Helium-4, Helium-3 and ultra-cold atoms. It has been widely used 
in making high-precision devices and also in cooling various systems. There have been extensive 
experimental search for possible exciton superfluid (ESF) in semiconductor electron-hole bilayer 
(EHBL) systems below liquid Helium temperature. Exciton superfluid are meta-stable and will 
eventually decay through emitting photons. Here we find that the light emitted from the excitonic 
superfluid has unique and unusual features not shared by any other atomic or condensed matter 
systems. We show that the emitted photons along the direction perpendicular to the layer are in 
a coherent state with a single energy, those along all tilted directions are in a two modes squeezed 
state. We determine the two mode squeezing spectra, the angle resolved photon spectrum, the 
line shapes of both the momentum distribution curve (MDC) and the energy distribution curve 
(EDC). By studing the two photon correlation functions, we find there are photon bunching, the 
photo-count statistics is super-Poissonian. We also stress the important difference between the 
quasi-particle excitations in an equilibrium superfluid and those in a stationary state superfluid. 
This difference leads to the explanation of recent experimental observation of excitation spectrum 
of exciton-polariton inside a planar cavity. We discuss how several important parameters such as the 
chemical potential, the exciton decay rate, the quasiparticle energy spectrum and the dipole-dipole 
interaction strength between the excitons in our theory can be extracted from the experimental 
data and comment on available experimental data on both EDC and MDC. We suggest that all 
the predictions achieved in this paper can be measured by possible future angle resolved power 
spectrum, phase sensitive homodyne measurements, and HanburyBrown-Twiss type of experiments. 
We demonstrate explicitly that the photoluminescence from the exciton in EHBL systems is a very 
natural, feasible and unambiguous internal probe of the nature of quantum phases of excitons in 
EHBL such as the ground state and the quasi-particle excitations above the ground state. These 
remarkable features of the photoluminescence can be used for high precision measurements, quantum 
communication, quantum information processing and also for the development of a new generation 
of powerful opto-electronic devices. 



I. INTRODUCTION 



An Exciton is a bound state of an electron and a 
hole. Exciton condensate was first proposed more than 
3 decades ago as a possible ordered state in solids 1 . 
Keldysh and Kozlov- 2 . argued that in a bulk semiconduc- 
tor, in the dilute limit n ex o? ex <C 1 where n ex is the ex- 
citon density and a ex is the exciton radius, the excitons 
behave as weakly interacting bosons, the exciton effective 
mass M is even smaller than an electron mass, for exper- 
imentally accessible exciton densities, the 3 dimensional 
Bose-Einstein condensation (BEC) critical temperature 
can be estimated to be ~ K. So in principle, the exci- 
tons can undergo BEC and become an excitonic super- 
fluid state below a few K. In the dense limit n ex af x ^> 1, 
the fcrmionic nature of the electrons and holes in the ex- 
citon will show up, the strong pairing BEC superfluid will 
crossover to weak pairing BCS superfluid—. However, in 
reality, it is very difficult to realize the BEC of excitons 
experimentally in a bulk system, because the short life- 
time r ex and the long lattice relaxation time tt, which is 
needed for the hot exciton gas to reach the cold temper- 
ature of the underlying lattice by emitting longitudinal 
acoustic phonons. Although exciton gas, bi-excitons and 
electron- hole plasma (EHP) phase have been observed 
in different bulk semi-conductors, no exciton superfluid 



phase has been observed in any bulk semi-conductors. 

Recently, degenerate exciton systems have been pro- 
duced by different experimental groups with two dif- 
ferent methods in quasi-two-dimensional semiconductor 
GaAs/AlGaAs coupled quantum wells structure^^^. 
When the distance between the two quantum wells is 
sufficiently small, an electron in one well and a hole in 
the other well could pair to form an exciton which be- 
haves as a boson in dilute exciton density limit ( Fig. la). 
This kind of inter-layer excitons are called in-direct ex- 
citons. In Butov and Snoke's labs ( also Bell lab 
the excitons are created by optical pumping and then a 
electric field is applied along the z direction to separate 
electrons from holes by a distance d ~ 30nm. There are 
also the current efforts from Bell lab£ which focused on 
the effects of electrostatic traps to confine the excitons 
in a given regime. In the undoped electron-hole bilayer 
(EHBL) sample prepared in Mike Lilly' lab^OLiid 2 . 
and the Cambridge group 13 ' 14 which is a heterostructures 
insulated-gate field effect transistors, separate gates can 
be connected to electron layer and hole layer, so the den- 
sities of electron and holes can be tuned independently 
by varying the gate voltages. Low densities and high 
mobilities regimes for both electrons and holes can be 
reached. Transport properties such as Coulomb drag can 
be performed in this experimental set-up. 
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FIG. 1: (a) In the experiments i n 4 ' 5 ' 6 , A laser with excita- 
tion power P ex is used to excite electrons from the valence 
band ( V ) to the conduction band (C) in the semiconduc- 
tor GaAs/AlGaAs electron-hole bilayer systems. Then a gate 
voltage V g is applied along the z direction to separate the 
electrons in the conduction band from the holes in the va- 
lence band. One electron in one quantum well and one hole 
in the other quantum well are bound to form an indirect ex- 
citon. In the experiments the excitons are generated by 

gate voltages, (b) The geometry of the photon emission from 
the indirect excitons 



The quantum degeneracy temperature of a two dimen- 
sional excitonic superfluid (ESF) can be estimated to be 
T5 X ~ 3K for exciton density n ~ 10 10 cto~ 2 and ef- 
fective exciton mass m = 0.22too where mo is the bare 
mass of an electron, so it can be reached easily by He 
refrigerator. It was established that the indirect excitons 
in EHBL has at least the following advantages over the 
excitons in the bulk: (1) Due to the space separation 
of electrons and holes, the lifetime T ex of the excitons 
is made to be 10 3 ~ 10 5 longer than that of those in 
bulk semi-conductors, now it can be made as long as 
microseconds. Due to the relaxation of the momentum 
conservation along the z direction, the thermal lattice re- 
laxation time T£ of the indirect excitons can be made as 
10 -3 that of bulk excitons, now it can be made as short 
as nanoseconds. So r ex 3> Tt is well satisfied even for 
the direct semiconductor such as GaAs. (2) Because all 
the electric dipoles are aligned normal to the 2d plane, 
the repulsive dipole-dipole interaction is crucial to stabi- 
lize the excitonic superfluid against the competing phases 
such as bi-exciton formation and electron-hole plasma 
(EHP) phase. So EHBL is a very promising system to 
observe BEC of in-direct excitons. There are two impor- 
tant dimensionless parameters in the EHBL. One is the 
dimensionless distance 7 = d/as ( &s is the Bohr radius 
) between the two layers. Another is r s where r s as is the 
typical interparticle distance in a single layer. Recently, 
one of the authors proposed that the EHBL maybe a 
more favorable system to observe a metastable excitonic 
supersolid (ESS) than the Helium 4 system^. The global 
phase diagram at T = labeled by the two parameters 
is shown in Fig. 3. 

Note that in all the previous experiments 4 *^ 6 -, a laser 
beam was consistently shined on a given 'bright' ring ( 
or a 'bright' spot ), however, the excitons will move to 
different locations which is at the center of the ring ( or 
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FIG. 2: A laser ( black arrow ) with excitation power P ex is 
used to excite electrons which are relaxed ( red arrow ) to the 
bottom of the trap, then emit photons ( blue dashed arrow ) 
spontaneously. 



a ring ) before they annihilate and emit lights ( Fig. 2 
). So in the stationary process of emitting lights, the 
number of the exciton condensation N is kept to be a 
constant. The laser beam was used to photo-generated 
the excitons, so it plays the role of a pump, however, 
because the exciton condensation happens at different 
locations than the laser pumping point, so if the BEC of 
excitons indeed happens near the center of the trap, it is 
indeed spontaneous instead of being stimulated. Indeed, 
as temperature^ is decreased from ~ 20K to ~ 1.7K, the 
spatially and spectrally resolved PL peak density center- 
ing around the gap in Fig. la E g ~ 1.545eV increases, the 
exciton cloud size decreases to L ~ 30/im, the peak width 
shrinking to ~ lmeV at the lowest temperature ~ 1.7K. 
All these facts indicate the possible formation of exciton 
condensate around 1.7 K. In our theoretical analysis in 
this paper, we assume the exciton cloud already reached 
the lattice temperature by interacting with lattice acous- 
tic phonons within the thermal lattice relaxation time T£ 
during its relaxation process to the bottom of the trap, 
it also become a superfluid through mutual dipole-dipole 
repulsive interaction and start to radiate photons at the 
exciton lifetime r ex tl, then we will calculate all the 
characteristics of the photons emitted from the exciton 
superfluid. We will compare our theoretical results with 
the experimental data in Sec. IX. The transient photo- 
luminescence from the excitons created by a short laser 
pulse will be discussed in a separate publication. 

In parallel to search for exciton superfluid in EHBL, 
extensive activitie s 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 ' 24 have also been 
lavished on searching for exciton superfluid in electron- 
electron bilayer system in the same semi-conductor ma- 
terial AlGaAs/GaAs subject to a high magnetic field 
in quantum Hall regime at total filling factor vt = \. 
When the interlayer separation d is sufficiently large, the 
bilayer system decouples into two separate compressible 
v = 1/2 layers. However, when d is smaller than a critical 
distance d c , the system may undergo a quantum phase 
transition into a novel spontaneous interlayer coherent 
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exciton superfluid phased. The exciton in this system 
can be considered as the pairing of an electron in top 
layer and the hole in the bottom layer after making a 
particle-hole transformation in the bottom layer. Other 
phases such as pseudo-spin density wave phase in some 
intermediate distance regimes was also propose d 22 ' 23 ' 24 . 
At low temperature, with extremely small interlayer tun- 
neling amplitude, Spielman et al discovered a very pro- 
nounced narrow zero bias peak in this possible exciton 
superfluid stated. M. Kellogg et al also observed quan- 
tized Hall drag resistance at h/e 217 . In the counterflow 
experiments, it was found that both linear longitudinal 
and Hall resistances take activated forms and vanish only 
in the zero temperature limit 1 ^. However, despite the in- 
tensive theoretical research 1 ^ in the past, there are still 
many serious discrepancies between theory and the ex- 
periments. It remains unclear if the excitonic superfluid 
was indeed realized in the BLQH system. 

It is instructive to compare the measurements to de- 
tect possible exciton superfluid in the BLQH and EHBL 
In the BLQH, there are mainly three kinds of transport 
experimental measurements (I) Quantum Hall resistance 
(2) Interlayer tunneling (3) Counterflow. In contrast to 
these quantum phases in BLQH which are stable ones, 
all these excitonic phases in EHBL in the Fig. 3 are just 
meta-stable states which will eventually decay by emit- 
ting lights. So the most natural experimental measure- 
ment for photo generated EHBL is the photolumines- 
cence (PL) which is quite different from all the transport 
measurements in the BLQH. The geometry of the photo- 
luminescence from EHBL systems is shown in Fig. Ib. In 
fact, the photon emission in EHBL plays a similar role 
as the interlayer tunneling the BLQH, so the theoreti- 
cal results achieved in both systems should shed on and 
transfer lights to each other. Very recently, transport ex- 
periments such as Coulomb drag were also performed in 
EHBL generated by gate voltagesZ&^ i 10 ' 11 ' 12 . It is pos- 
sible to also perform counterflow experiment in the near 
future. The PL experiment can also be performed in 
this gate voltage generated EHBL, although the emitted 
lights are weaker than those from the photo generated 
EHBLS. 

In parallel to the experimental search for the ex- 
citon superfluid in the EHBL, there are also ex- 
tensive experimental activities to search for exciton- 
polariton superfluid inside a planar micro-cavity. Al- 
though exciton condensation in a single quantum well 
( SQW ) has not been observed so far, there are 
some evidences for the observation of Exciton-polariton 
(EP) condensation in SWQ enclosed inside a planar 
microcavity^ 3 ^ 3 ^ 3 ^^^^ 2 .. These evidences 
include macroscopic occupation of the ground state, spec- 
tral and spatial narrowing, a peak at zero momentum in 
the momentum distribution ( see Fig. 8b ) and sponta- 
neous linear polarization of the light emission and so on. 
The elementary excitation spectrum of exciton-polariton 
was also found^i to be very similar to that in an equi- 
librium superfluid with notable exceptions near k = 0. 



d/a B 





Strong 

?££)\ < BEC » 

Excitonic supersfluid 

(b) 



FIG. 3: Zero temperature phase diagram of electron and hole 
bilayer at the balanced case n e — nh ( but m e 7^ mj, ). The 
very important effects of spins, traps and disorders are not 
taken into account. The d is interlayer distance. r s ag where 
as is the Bohr radius is the average inter-particle distance in 
a single layer. At high densities ( low r s ), along the dashed 
line (I), there is a transition from the weak pairing ( or BCS 
) ESF to the electron-hole plasma ( EHP ) phase. At low 
densities ( high r s ), along the dashed line (II), there is a 
transition from the strong pairing ( or BEC ) ESF to the 
excitonic supersolid (ESS), then to excitonic normal solid ( 
ENS ). Inside the ESF, there is a BCS ( weak-pairing ) to 
BEC ( strong pairing ) cross-over tuned by r 3 denoted by the 
dashed line. The EHP phase is a conducting phase, while the 
ENS is an insulating phase, so the EHP to ENS transition can 
be considered as a metal to insulator transition. Thick line is 
1st order transition, the thin one is second order. The EHP 
to the weak pairing BCS transition may also be 1st order. 



This puzzle will be resolved in Sec.VII-6. Recently, 
several ultra-cold atom experiment o 43 ' 44 ' 45 ' 46 successfully 
achieved the strong coupling of a BEC of N ~ 10 5 87 Rb 
atoms inside a cavity. Motivated by these achievements 
of SQW and atomic BEC embedded inside a cavity, we 
suggest that in near future experiments, the EHBL can 
also be enclosed in a planar micro-cavity, so one can 
search for possible superfluid of indirect exciton polar- 
tion (IEP). One advantage of the EHBL over the SQW 
is that as shown in the previous paragraph, the dipolcs 
of the indirect excitons are all aligned along the z direc- 
tion, so the dipole-dipole interaction is repulsive, this also 
guarantees the IEP-IEP interaction is repulsive which is 
a sufficient and necessary condition to stabilize a super- 
fluid against other possible states. This strong coupling 
regime in a planar micro-cavity will be investigated in 
separate publications. 

Although there exist extensive experiment measure- 
ments on photoluminescence from presumably achieved 
exciton BEC in the electron-hole bilayer (EHBL) 
system^i^, so far, there is no systematic theory on how 
photons interact with the indirect excitons in different 
quantum phases in the EHBL system and how the charac- 
teristics of photons can reflect the nature of the quantum 
phases in the Fig. 3. In this paper, we will study three 
dimensional photons interacting with the two dimension 
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indirect excitons in the excitonic superfluid phase in the 
BEC side in the EHBL. We will work out how the pho- 
toluminescence from this phase can reflect the properties 
of both the condensate and the Bogoliubov quasi-particle 
excitations above the condensate at zero temperature 
T = 0. We find that due to the non- vanishing order pa- 
rameter in the ESF phase, the emitted photons along the 
direction perpendicular to the layers ( namely with zero 
in-plane momentum ) are in a coherent state, while the 
non- vanishing anomalous Green function in the ESF lead 
to a two mode squeezed state of the emitted photon along 
all tilted directions ( namely, at finite in-plane momenta 
) as shown Fig. fb. We determine the angle resolved 
power spectrum, squeezing spectrum, one and two pho- 
ton correlation functions along all the possible directions 
including normal and tilted directions. From the two 
point correlation function, we can identify the quantum 
nature of the emitted photons such as photon bunch- 
ing, ant i- bunching, also the photo-count statistics such 
as super-Poissonian, Poissonian and sub-Poissonia n 47 ' 48 . 
We will also determine the momentum distribution curve 
(MDC) and energy distribution curve (EDC) which are 
the integrated angle resolved power spectrum at fixed 
energy and fixed momentum respectively and then com- 
pare with the available experimental data on EDC. We 
will suggest that all our predictions can be measured 
by possible future angle resolved power spectrum, phase 
sensitive homodyne measurements, and HanburyBrown- 
Twiss type of experiments. We will also elucidate the 
physical reasons why the angle resolved power spectrum 
takes the super-radiant form even in the thermodynamic 
limit when the exciton decay rate is sufficiently large, why 
the characters of the light emitted from the ESF phase 
can reflect both the nature of the ground state and the 
Bogoliubov quasi-particle excitations above the ground 
state even at T = 0. The photoluminescence from the 
other phases in EHBL system will be studied in subse- 
quent works. In this paper, we did not consider the very 
important effects of spins of electrons and holes which 
lead to the formation of the bright excitons with J = ±1 
and the dark excitons with J = ±3/2 26 , the effects of the 
trap, finite thickness and disorders in the sample^ 7 .. All 
these will be discussed in separate publications. 

The rest of the paper was organized as follows, in sec- 
tion II, we will derive the interaction between the 3 di- 
mensional photons with 2 dimensional indirect excitons 
in the exciton gas phase in EHBL with a distance d. In 
section III, we will derive the total Hamiltonian in the 
ESF phase on the BEC side, separating the interaction 
into the coupling to the condensate at zero in-plane mo- 
mentum k — and to the Bogolubov quasi-particle at 
non-zero in-plane momentum k ^ 0. Then we will show 
how a coherent state is emitted at k = and discuss 
several remarkable properties of power spectrum along 
the normal direction in section IV. In Section V, we will 
develop systematically an input- output formalism for a 
stationary state. Then by using the input-output formal- 
ism, we will calculate various emitted photon character- 



istics in the follow sections. In Section VI, by calculating 
the squeezing spectrum, we show that the ESF phase of 
the excitons play a similar role as a two mode squeez- 
ing operator which squeeze the input vacuum state into 
a two mode squeezed state, so the emitted photons at 
non-zero k are always in a two mode squeezed state even 
off the resonance. In section VII, we evaluate the angle 
resolved power spectrum, the line shapes of both MDC 
and EDC. We find that the angle resolved power spec- 
trum takes a stationary super-radiant form even in the 
thermodynamic limit when the exciton decay rate is suf- 
ficiently large compared to the energy of the Bogoliubov 
excitation. We also compare with the Dicke model on 
super-radiance of finite N atoms in conventional quan- 
tum optics. By working out the special nature of ex- 
citation spectrum in an non-equilibrium superfluid, we 
resolve the puzzle observed in^ 7 -. We also resolve the In 
section VIII, we compute the one and two photon corre- 
lation functions at non-zero in-plane momentum k 
and find the photon statistics at any non-zero k. In sec- 
tion IX, we will compare our theoretical results on EDC 
and MDC achieved in the last section with the previ- 
ous experimental PL data in 4 -^ and also discuss possible 
future experimental set-up such as angle resolved power 
spectrum measurement, phase sensitive homodyne mea- 
surement, and HanburyBrown-Twiss type of experiments 
to test the predictions achieved in sections IV- VIII. In the 
final section X, we summarize the main results on the co- 
herent, squeezed and macroscopic super-radiant nature of 
the emitted photons and point out their crucial difference 
than the previous coherent and squeezed states generated 
by pumps in non-linear media. In the appendix A, we will 
explain how the exciton superfluid emit photons in terms 
of a intuitive Radiation Zone picture. In the appendix 
B which supplement section VI, we give a more intuitive 
proof that all the emitted photons along the tilted di- 
rections are in a two mode squeezed state even off the 
resonance. In the appendix C, we clarify the relation 
between the quantities calculated in the main text and 
experimental measurable quantities. In the appendix D, 
we will perform a Golden rule calculations to second or- 
der at both k = and k ^ by using the many body 
exciton BEC ground and excited states with Bogoliubov 
quasi-particles and compare with the results achieved in 
section IV by Heisenberg equation of motion at k = and 
in section VI- VIII by the non-perturbation input-output 
formalism calculations at k =/= 0. 

II. THE INTERACTION OF EXCITON WITH 
PHOTON IN THE BEC SIDE OF EHBL: A 
MICROSCOPIC POINT OF VIEW 

In this section, we will derive the coupling constant 
between the three dimensional photons with two dimen- 
sional indirect excitons from microscopic point of view. 
The second quantization Hamiltonian consists of three 
parts H = Ha + H ex + Hi nt where the first part is the 



5 



Hamiltonian of free photons: 



Ha — y^^fcQ fc-A afc,A 

k,X 



(1) 



e k (^fc) * s annihilation operator of the electron (hole), 
then the cxciton Hamiltonian can be rewritten as 



E v {k)hih z ] 



(5) 



where a^A (^1 A ) is the annihilation (creation) operator 
of the photon, it has polarization A and three dimension 
momentum k = (fc, k z ) where k is the two dimensional in- 
plane momentum. The frequency of the photon is LOk = 

Vgtjk 2 + k 2 , where v g — c/^/e. Here, c is the light speed 
in the vacuum and e ~ 12 is dielectric constant of GaAs. 

In EHBL, we can decompose the electron field into two 
parts: 



(2) 



where f stands for two dimensional positions in the 
EHBL, $i(z) and $2(2) are strongly localized around 
z = and z = d respectively. Then ipi(r, z) ~ z — 
0) = ipi(r) and ip2(r, z) ~ ^(r, z = d) = 4>2{r) are elec- 
tron operators in top and bottom layers in Fig. la. The 
second part is the Hamiltonian of the exciton: 



Hq + V int 



d 2 r4{r)[ 



2rriQ 
7i 2 V 2 



2mo 



V in t = -JdWSpmVi^-r^ir') (3) 

where V\ c and V 2c are periodic crystal potentials in the 
two layers, 6pi(r) = ipi(r)ipi(r) — rii,i — 1,2 are normal 
ordered electron densities on each layer. The intralayer 
interactions are Vii = V22 = e 2 /e|r|, while interlayer 



interaction is V12 = V21 = e 2 /ey\r\" + d 2 where e is the 
dielectric constant. 

Considering the effects of the crystal potentials V\ c and 
Vic in the two layers, the electron field operators in the 
two layers can be expanded in terms of Bloch waves: 

fort = -^Y, u c,k^zi = oyZ%^ 

k 

= j= E = dV k -% h (4) 

where S is the area of the layers and the Bloch wave func- 
tions satisfy H^u i ^(r)e ik ' r — EiU i k "{r)e lk r with i = c, v. 
For direct semi-conductor such as GaAs, there is a min- 
imum ti in the conduction band i = 1 and a maximum 
e v at the valance band i — 2. We will set e v = below, 
then the band gap is E g — €4 — e v — e v . 

It is convenient to perform the particle-hole transfor- 
mation in the valence band c. ? = er and c„r = h} -., then 

lk k 2k 



In the rest of the paper, we consider the dilute limit 
along the path II in Fig. 3 where the size of the exciton is 
of the order of distance between the two layers. The in- 
teraction between electron and electron (or hole and hole) 
in the same layer will just renormalize the masses of the 
electron (hole) in the same layer—. Then the Hamilto- 
nian of the exciton can be further simplified to 



H R 



= E[(^r + ^h+^r4^ + K-„ (6) 



£ k 2m/j * k 



where m e (m/j) is the effective mass of the electron 
(hole). The exciton creation operator is defined as 



bl = J^p^oip — rriek/M)e~hX where the exciton mass 

M = m e + rrih and fo(p) is the Fourier transformation 
of the wave function <po(r) which satisfies 

,,2 



1 



2m. 



-E b cj> (r) (7) 



eJ\r\ z + d 2 



where the reduced mass is l/m r = l/m e + l/rrih- 

For a direct exciton in a single quantum well, d = 
in Eqn[7l the size of an exciton is a ex — h 2 e/e 2 m r — 
e^ffls ~ iOOas ~ 50A where a B = h 2 /e 2 m ~ 0.53A 
is the bare Bohr radius with the S'-wavefunction 4>o {f) = 
(^2-) 1 / 2 e~ 2r / a " x . The binding energy was known to be 



E b 



-e 2 /2a e 



-lOmeV. For an in- 



niQ e 2 2o,b 

direct exciton in the EHBL, d > 0, the exact form of 
the solution of Eqn[7] is not known, but it is not needed 
in the following discussions. Taking the exciton density 
n ~ 10 10 cm~ 2 , we can see that the average spacing be- 
tween excitons a ~ lOOnm 3> a ex ~ 5nm, so the sample 
is in the dilute limit. The exciton operators satisfy the 
commutation relation [b%, bV\ = 1 approximately in the 
dilute limit along the path II in Fig. 3. This approxi- 
mation is valid when the electron and hole form a tight 
bound state, the pair breaking process into electron and 
hole is at very high energy and can be neglected at low 
temperature, the excitonic system is essentially a bosonic 
system. Finally the Hamiltonian of the free exciton reads 



(8) 



where Ef = h z k 2 /2M + E g - E b . The third part is the 
interaction between excitons and photons which can be 
separated into one photon and two photon parts: 



H 



(a) 



e/mo / d 3 r^l(r)A(r) • pip(r), 



H 



= e 2 /2m / dTripT(r)A 2 (r)ip(r), 



(9) 
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where r = (r, z) stands for the three dimensional position, 
mo is the bare mass of an electron, the vector potential 
of the photon is: 



k ■ r 



t -ik-r 
l k,X e 



) (io) 



A'. A 



where V = L 2 x L z is the normalization volume of the 
whole 3-dimensional system and the Zk\ is polarization 
of the photon with three dimensional momentum k. 
By inserting the vector potential A(r) and the electron 

(hole) field into the interaction Hamiltonian H^ t , the 
approximate relation^ J d 2 fu* (f)u Vt o(r) ~ leads to 

H$ t ~ in the Hilbert space of the exciton. Then we 

can project the interaction Hamiltonian H^ t into the 

Hilbert space of the excitons 



6t|0. 



where the 



k is the two dimensional momentum of the exciton. In 
this subspace of the exciton, the interaction Hamiltonian 
is 



H 



(a) 



Ei 

k 

Ei 

k 



Uk ex H^ t \0 ex ) k ex ){0 ex \+h.c] 



Hg> |0e 



(11) 



By utilizing the electron field operator and the re- 
lation [r,Ho] = ip/riiQ, we find the matrix element 

^int |0«b) = *Sfc ff(^)°fc where the coupling con- 



kex 

stant is 

g(k) 



-{l/2eu k L z ) l l 2 Efu. c 



,{k)m~L- z 1/2 as) 



where L z is the normalization length along the z direc- 
tion, Hcv(k) — ek\ ■ Dk where the transition dipole mo- 
ment between the conduction band and the valence band 
is D k = (D xy ,D z ): 

D k z = Jd 3 re ik 'X/r,z)ezu Vi£ (f,z)^ 1 (z)^(zlU) 

Both are essentially the overlap between the wavefunc- 
tion of the electron in conduction band in one quantum 
well and the wavefunction of the hole in valence band in 
the other quantum well which lead to a small interlayer 
tunneling. 

For a given photon momentum with k, the polarization 
e*fe2 = k x Dk in Fig. lb is normal to the transition dipole 
moment, so can be dropped out, we need only consider 
the single polarization tki which is in the plane deter- 
mined by k and Dk in Fig. lb, then fx cv (k) = Dksm8k- 
Note that the transition dipole moment Dk from the 
conduction band to the valence band at a momentum 
k is completely different from the static dipole moment 
D s = e J d 2 r\\ip2\ 2 — iV'il 2 ] i n the dipole-dipole interac- 
tion Vd(q) in Eqn. [15j Although D s is completely along 



the z-direction in the dilute limit along the path II in 
the Fig. 3, the Dk is along a general direction depending 
on k shown in Fig. lb. For example, in the absence of 
interlayer tunneling, Dk = 0, but D s ^ 0. 

Finally, the interaction Hamiltonian is simplified to 



H^ t =J2Mk)akbl + h.c.}, 



(14) 



where there is a in-plane momentum comservation be- 
tween emitted photons and the excitons. 

In the summary of this section, we derived the Hamil- 
tonian H = Ha + H ex + Hfol of the indirect exciton 
+ photon Hamiltonian in the EHBL system given by 
Eqns H18114I 



III. THE COUPLING BETWEEN THE 
PHOTON AND THE CONDENSATE, THE 

PHOTON AND THE BOGOLIUBOV 
QUASI-PARTICLES IN THE EXCITONIC 
SUPERFLUID 

In this section, we will consider the effective interaction 
between the photon and the Bogoliubov quasi-particle 
excitations in the excitonic supcrfluid phase in the BEC 
side in Fig. 4. The total Hamiltonian in grand canonical 
ensemble is the sum of excitonic superfluid part, photon 
part and the coupling between the two parts H t = H s f + 
H ph + H int where : 



Hsf 
Hph 

Hint 



k kpq 

k 

J2Mk)a k bt + h.c], (15) 



where Ef = h z k 2 /2M + E g — E b , S is the area of the 
sample, V d (q) 



1 — e q ) is the dipole-dipole inter- 
action between the excitons^, Vd(|f| ^> d) — e 2 d 2 / \ f\ 3 
and Vd(0) = 2ne which is a finite constant leading to 
a capacitive term for the density fluctuation^. It is im- 
portant to stress that in a stationary state, the chemi- 
cal potential p for the excitons in Eqn 11 51 is kept fixed 
by the off-resonant pumping which is the laser pump- 
ing iniiSiS 7 and the gate voltage pumping i n 7 ' 8 ' 9 ' 10 ^ 11 ' 12 . 
Very similar point was also stressed in^ in the context of 
non-equilibrium stationary transport through a quantum 
dot. 

In the dilute limit, Vd is relatively weak, so we can ap- 
ply standard Bogoloubov approximation to this system, 
this is in contrast to Helium 4 system which is a strongly 
interacting system. The phase 6 representation used ir>i£ 
is very useful to study vortex anti- vortex excitations and 
Kosterlitz-Thouless transition at finite temperature. In 
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this paper, we focus only at T = 0, so we can ignore 
the topological excitations in the phase winding 9 and 
just use the Bogoloubov approximation to treat the non- 
topological low energy excitations. So in the ESF phase, 
one can decompose the exciton operator into the conden- 
sation part and the quantum fluctuation part above the 
condensation 6^ = \fN8t Q + b^. Note that we ignored 
the zero point fluctuation above the condensate which 
is justified in the thermodynamic limit. In the real ex- 
perimental situation where finite number of excitons are 
trapped inside a trap, its importations will be addressed 
in a separate publication^. 

In the stationary state, the chemical potential is fixed 

at 



H = + nV d {0) = (E g - Et) + nV d {0) 



(16) 



which is determined by eliminating the linear term of 
bo in the Hamiltonian H s f. So the chemical potential 
is increased ( or called " blue shifted" ) from the single 
exciton energy E g — Eb by the dipole-dipole interaction 
nVd(0). In the experimental set-up shown in the Fig. la, 
E g - E b = E° g - E b - eEd where E° - E b is the bare 
conduction-valence band gap at zero gate voltage V g = 
and E is the electric field due to the applied voltage V g in 
Fig. lb. The exciton density n is determined by the laser 
excitation power P ex . So the chemical potential can be 
tuned by the two experimental parameters V g and P ex . 

Then the Hamiltonian of exciton BEC upto the 
quadratic terms is 

Hsf = 5> S + VS)n)l^b % + (^^SJSLs + h.c.)}, 

k 

(17) 

where the density of the condensate n = N/ S. For study- 
ing the quasi-particle excitation spectrum of the exciton 
BEC, we utilize the Bogoliubov transformation 



(18) 



to diagonize the Hamiltonian H s f, where the transfor- 
mation coefficients arc 



in terms of the quasi-particle creation and annihilation 
operators 3^ and 0l and E(0) is the condensation energy. 
The spectrum of quasi-particle excitation is 



E(k) 



2nV d (k)\, 



(21) 



As k -> 0, E(k) 
particle is: 



u k where the velocity of the quasi- 



u = y/nV d (0)/M 



2ne 2 dn 
eM 



(22) 



Plugging n ~ 10 10 cto~ 2 , d ~ 30nm,M ~ 0.22mo into 
Eqn l221 we find u ~ 5 x 10 cm/s. The quasi-particle 
spectrum is shown in Fig. 3 where the roton mode is 
due to the long-range dipole-dipole interaction^. Even 
at T = 0, the number of excitons out of the condensate 



d 2 k 2 
(2tt) 2 ^ 



(23) 



which is the quantum depletion of the condensate due to 
the dipole-dipole interaction. From Eqnfini we can see 
w| — > 1/k, as k — > and — > 1/fc 6 , as k — > oo, so Eqnl23l 
is well defined. 

For studying the effects of the condensate and the 
quasi-particle excitation to the emitted light separately, 
we decompose the interaction Hamiltonian Hi into the 
coupling to the condensate part 



Na 



k, + 



h.c], 



(24) 



and the coupling to the quasi-particle part: 
H ! = J2i i 9(k)a k bt + h.c.}. 



(25) 



nV d {k) 



2E(k) 
+ nVd(fc) _ 1 
2E(k) 2 



1 

2' 



1 and u^v^ 



„ which is completely 

k k k k 2E(k) 1 J 

due to the exciton dipole-dipole interaction. The quasi- 
particle creation and annihilation operators 3r and 31 

fc k' 

satisfy the Bose commutation relation [/%,/?£,] = 
Finally, the Hamiltonian of exciton BEC is given by 



In the following two sections, we discuss the properties of 
emitted photons with the zero in-plane momentum k = 
and the non-zero in-plane momentum k 7^ respectively. 
mq\ Note that due to the electron-hole asymmetry, we do not 
expect there is a up-down Z2 symmetry. However, for 
the simplicity of notations, we assume there is such a Z2 
symmetry, so we can treat the radiations in the upper and 
down half space on the equal footing. All our calculations 
can be generalized straightforwardly to take into account 
the asymmetry quantitatively in the real EHBL system. 

In the following, in order to keep the relative energy 
difference between the exciton and the photon intact, we 
made a rotation a — a(t)e~ v *, so we can focus on the 
(20) slowly varying d(t) and neglect the~ in the following sec- 
tions. 
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E(k), 




Linear mode 



FIG. 4: The spectrum of the quasi-particle in exciton BEC 
where there exist the linear mode and roton mode. This spec- 
trum only holds when teh lifetime of the exciton t £X — > oo. 
The spectrum will be changed to Fig. 6b when r ex is large, but 
finite. 



IV. THE COHERENT STATE, LINE WIDTH 
AND POWER SPECTRUM IN THE NORMAL 
DIRECTION AT k = 0. 

Because the condensate carries no in-plane momentum 
k = 0, so from H? nt in EqnlMl the Heisenberg equation 
of motion of the photon annihilation operator is 



id t a kz = (iJk z - M - i^Ofc. - ig*(k z )VN + iF(t) (26) 

where we have dropped the zero mode fluctuation of the 
zero momentum condensate which is negligible in the 
thermodynamic limit. The k is the decay rate of the pho- 
ton due to its coupling to a reservoir ( or bath ) , its value 
and physical meaning will be determined self-consistently 
in the following. The F(t) is the fluctuation of the reser- 
voir satisfying (b\F(t)\b) = where (0) b = Tr(p b O) 
denotes the mean value of the operator O at the bath 
state and pb is the density matrix of bath. From Eqnl2"6l 
one can see the exciton condensation y/N plays the role 
of an effective pump on the photon part. As shown in 
the Eqn|H g(k z ) ~ L~ 1/2 . 

Following the standard laser theory, we decompose the 
operator as its mean plus its fluctuation: a kz = ( a k,) in + 
Sa kz where the initial state |). in Fig. 2 is taken to be 
\BEC) \0) ph \b). Here, \BEC) denotes the ground state 
of the condensation; the initial photon state |0) h satisfies 
a kz (0) \0) p h = where a kz (0) is the photon operator a kz 
at the initial time t — 0; \b) denotes the state of the 
photon reservoir. From Eqn l261 it is easy to see that 

idt (a kz ) = (uj kz - [i- in) (a kz ) - ig*(k z )VN. 
The stationary solution is 

ig*(k z )VN 



(27) 



(uo kz — fi — m) 



(28) 



which is the photon condensation induced by the exciton 
condensation at k = 0. So the output state along the 
normal direction is a coherent state. 
The fluctuation Sa kz obeys: 

id t 5a kz = (uj kz - fi - in)5a kz + iF(t), (29) 

When t >> 1/k, the photon fluctuation is completely 
determined by the fluctuations of the reservoir: 



6ai 



Jo 



'k z —n){t—T)—n{t—T) 



(30) 



Where the fluctuation-dissipation relation dictates that 
(FHt)F(t>)) - Kn^J(t-t'),(F(t)FHt')) = n(n^ + 
l)6(t-t'),(F(t)F(t')) = (F^t)F^t')) = 0. 
From Eqn l301 we find: 

(5al z (t)6a kz (t')} = ^nk^-^e^ -/*)(*-t') (31) 

The power spectrum of the fluctuation is given by the 
Fourier transformation of the correlation function: 



S{w) 



Sal z (t)5a kz (t')\ e-^-^ t - t ">dt (32) 



where we have defined the photon frequency with respect 
to the chemical potential fi. 

By inserting Eqn. EUinto EqnJ351 we get and the power 
spectrum 



1 



2(w-w fe J 2 + K 2 ' 



(33) 



where the particle distribution of the photon reservoir is 
n^j = l/(e^/ T — 1) and T is the temperature of the reser- 
voir. The result Eqnj33] is consistent with the Wiener- 
Khintchine theorem. 

The number of the emitted photon is 



V = \ a k z a kz ) = 



N\g(y/c)\ 2 

(Ul kz - /i) 2 + K- 



+ ^/2, (34) 



where we have set g*(k z ) around uj kr = p. Because the 
condensate N is much larger than the particle distribu- 
tion of the bath n UJhz , so when the temperature of the 
reservoir T — * which is the case considered in this pa- 
per, the number of the emitted photon is dominated by 
the first term. 

Now we will determine the value of k self-consistently. 
The total number of photons is 



iV = 5>^ =N{\g\ 2 D)/n 



(35) 



where D — L z /v g is the photon density of states at k = 0. 
Note that the exciton decay rate 70 = \g\ 2 F> at k = is 
independent of L z ( see Eqn. |44] and [48] for general ex- 
pressions of the density of state D^{oj k ) and the exciton 



9 



decay rate 7^ at any k ). The N p h has to be proportional 
to L z in order to get a finite photon density in a given vol- 
ume L 2 x L z in a stationary state. This self-consistency 
condition sets k = v g / L z — > so that 

N ph = N lQ L z /v g ~ L z (36) 

Plugging this value of k into EqnJM] leads to: 



(37) 



which is independent of L z as required ! We showed that 
the power spectrum emitted from the exciton condensate 
has zero width. This conclusion is robust and is indepen- 
dent of any macroscopic details such as how photons are 
coupled to reservoirs. 

The radiation rate from the condensate is: 



PZ d = iv 7oM 



(38) 



V. THE INPUT-OUTPUT FORMALISM FOR A 
STATIONARY STATE AT k / 

In this section, we consider the photons with in-plane 
momentum k ^ 0. From Eqnf25l it is easy to see that due 
to the in-plane momentum conservation, the exciton with 
a fixed in-plane momentum k coupled to 3 dimensional 
photons with the same k, but with different momenta k z 
along the z-direction, so we can view these photon acting 
as the bath of the exciton by defining = ^ fc g{k)a k . 
As shown in Eqn l23l due to the dipole-dipole repulsion, 
even at T = 0, there are also excitons depleted from the 
condensate. These excitons will emit photons at non- 
zero k. By using the standard input-output formalism 
for a stationary state discussed in^i, we will achieve the 



squeezed spectrum, angle resolved power spectrum and 
photon correlation functions of the emitted photon in the 
following sections. 

The Heisenberg equations of motions of the photons 
and excitons arc 



d t a k = -i{uj k - fi)a k - g*(k)b£, 
d t B % = -iZB % + A % , 

where B % = (l^f, A % = (1^1* )* and 



(39) 



Ti= ( e % + nV d {k) nV d (k) \ 

\ -nV d (k) -e % -nV d {k))' V 1 

The formal solution of a k can be written either as the 
initial state at to < t or the final state at t\ > t: 



which is also independent of L z as required. In fact, in 

the limit k = v g /L z — > 0, Eqns. [33] becomes S(u>) = 

n Uk S(u> — tu kz ) which is negligible at low temperature. and 

From the Eqnf571 one can see that both the pumping 

term ig*y/N ~ L z *' 2 and the dumping term in ~ L^ 1 

approach zero as L z — > 00 limit in such a way that a 

stationary state is reached. 

From the experimental data in section IX, taking N ~ 
10 5 ,7o ~ QA(j,eV,fi~ 1.54eV, we find P^ d ~ lfiW. 

In summary, the coherent light emitted from the con- 
densate has the following remarkable properties: (1) 
highly directional: along the normal direction (1) highly 
monochromatic: pinned at a single energy given by the 
chemical potential fi (3) high power: proportional to the 
total number of excitons. These remarkable properties 
are independent of any microscopic details as such as the 
excitation power P ex and the line width of the pump- 
ing laser as long as they can generate excitons across the 
band gap. This fact could be useful to build highly pow- 
erful opto-elcctronic device. In the appendix D, we will 
give a more intuitive derivation of these results from a 
golden rule calculation. 



g*(k) / dt'lAt^e-^-^-^ (41) 



+g*(k) / dt'b^e-^*-^-^, (42) 



When plugging Eqns. |4"T1 and |4"21 into EqnHM we find 
it is convenient to define the input and output fields as: 

af(t) = V ; 1 a k (t Q )e^ {uJk -^ {t ~ f0 \ 

°r(*) = ~E ,nt , °*(*i)e- i( "*- /,)(t - tl) ,(43) 

where the density of states of the photon with a given 
in-plane momentum k is 



D k {u k ) 



UkL z 



V 9\l"t-V 2 g 



(44) 



which is proportional to L z . It can be shown that if 
v g \k\ <C H, the input and output fields obey the Bose 
commutation relations: 

[4"(t),a£t(t')] = [af^a^it 1 )] = 5, p 5(t-t') (45) 



In term of the input field sf(t) = {af(t), a m l(i)) T , 
the exciton operator B k obeys 



9 t fl s =(-ffi-|) J B £ + V7 £a|', 
In terms of the output field a°^ lt (t) 



(46) 



(a|«'(i), a°J{t)) T ,it obeys 



' >k A k ' 



(47) 
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where the effects of photon-exciton coupling are com- 
pletely encoded in the exciton decay rate 7^ = 



7fc 



Ef 2 \D k \\m 2 e k tf(0) 



U/2 



(48) 



where 9 k is the angle between the transition dipole mo- 
ment D k and the 3 dimensional photon vector k shown in 
Fig. lb. Note that 7^ is independent L z , so is an experi- 
mentally measurable quantity as shown in section IX- 1. 
From the rotational invariance in the Fig. lb, we can con- 
clude that 7£ ~ const. + \k\ 2 as k — > as shown in Fig. 6.. 

When comparing Eqn l46l and 1471 with Eqn l26l we can 
see that input photon field ^/T^al™ ( or the output pho- 
ton field ^/T^ajp* ) which are summation of continuous 

spectral of photons at a given in-plane momentum k, but 
with different k z as shown in Eqnf43] plays the role of the 
reservoir F(t) for the excitons, while the exciton decay 
rate 7g due to the photon-exciton coupling plays the role 
of k. However, there is no similar source ( or pumping ) 
term like — ig*{k z )\fN . 

When deriving EqnHB] and 271 we have assumed the 
density of state D^(uj k ) at a given in-plane momentum k 
varies slowly around the characteristic frequency u>k = A*- 
Indeed as shown later in Figs. 7. 8. 10, the maximum of 
squeezing spectrum and power spectrum is very narrowly 
peaked around to k = fx, so it is reasonable to set uj k = \i 
in DAwk)- This approximation is essentially a Markov 
approximation which is valid only when the in-plane mo- 
mentum k is much smaller than fi/v g in Eqn l48l In fact, 
as to be shown in section VI- 1, the maximum in-plane 
momentum k max = fi/v g . This is also the same ap- 
proximation for the commutation relations Eqnl45l hold. 
However, when the emitted photon is along the plane, 
namely with k z = 0, the Markov approximation becomes 
in- valid. So all our following calculations are valid as long 
as the emitted photons are not too close to along the xy 
plane. 

The relation 



in 



(49) 

is derived from Eqn l46l and Eqn l47l The Fourier trans- 
formations of Eq. (gU), Eq. (g7|) and Eq. (gSJ) lead to 
input-output relation: 



a<? t (cj) = W- 1 (u)W*(u)a™(Lj) 



where W(u>) 



-7gJ/2 + 



S) and a| ut (a/ 



(50) 
and 



aSy(u)) are the Fourier transformation of a- ,(t) and 

k v ' k x 1 

at ra (i). Then the component aS"*(o;) of the output field 
a J J,t (u) is related to the input fields by: 



af\u) = [-l + 7jg G n (fc> + i^)]af( W ) 



+7fcG a (fc,w + i— ) 



2 

Ik \ _mf 



where the normal Green function G n (k,u>) and the 
anomalous Green function G a (k,uj) arc 



G n (k,u>) = i 
G a {k,u) = - 



u> + eg + nV d (k) 
uu 2 -E 2 (k) ' 
inVd{k) 



E 2 (k)' 



(52) 



which are determined by the properties of the quasi- 
particle in the exciton BEC. In fact, they are just the 
retarded Green functions after making the analytic con- 
tinuation iui n — + u) + iS in the corresponding imaginary 
time Green functions. The exciton decay rate 7^ in the 
two Green functions in Eqn l51l just stand for the fact 
that the excitons are decaying into photons. Note that 
the Fourier transformation of the Eq. (|4"3"| leads to 



LU k - (J, 



(53) 



In fact, Eqnl51lcan be viewed as a S matrix relating the 
input photon field at to — > —00 to the output photon field 
at t\ — > 00. In the following sections, we will calculate 
the squeezed spectrum, angle resolved power spectrum 
and photon correlation functions of the emitted photons 
respectively. . 



VI. TWO MODE SQUEEZING SPECTRUM 
WITH k 

Eqn. I5T1 suggests that the output field is the two mode 
squeezed state between k and —k, so it is convenient to 
define A° t ut (uj) = [a° t ut {w) ± a ou Muj)}/y/2 and Af(uj) = 

[af{uj) ± a™Ju)]/V2- Then: 

A i%(u) = [-l + 7 % G n (k,u; + i^-)}AfJuj) 

±^ k G a (k,u + i^)Afl ± (-u>) (54) 

The position and momentum ( quadrature phase ) op- 
erators of the output field is defined by 

X± = A^Uuy^W + A^U-^e'^-^ 



fc,± 



iY± = A% u l(uj)e i ' t '±^ -^ u *t(_ w ) e -^±(-u.) ( 55 ) 



fc,± 



The squeezing spectral which measure the fluctuation 
of the canonical position and momentum are defined by 



% X± M = (X±(u)X ± {-w)) in , 

s UiY ju) = {y±(«)y±H) te . 



(56) 



where the 6(li + u>') function was omitted for notational 
simplicity, the in-state is the vacuum state of the in- 
put field \BEC)\0) shown in Fig. 4. Because the av- 
erage (X±(uj)} m = (Y±(w)) in = 0, then the squeez- 
(51) m g spectrums are Sx±(w) = \AX±(uj)\ 2 and Sy±(w) = 
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FIG. 5: The two mode squeezed state in the quadrature phase 
space. Only when the photon frequency resonates with the 
quasi-particle energy at E(k) > 7fc/2, the two mode squeezed 
states are mode + and mode -. All the other cases have non- 
zero squeezing angles shown as mode~±. See also Fig. 9. 



|AY±(a;)| 2 . It can be shown that 



l + {:X ± (u)X ± (-u) 



(57) 



where : AB : denotes the normal order of the A and B 
with respect to the A™*{iS), but the average is taken 

with the incoming vacuum state. 

For notational conveniences, we set <^_(u;) = 7r/2 + 
4>+{u)) and just set 4> + {uj) = <P{uj). Then we find 
Sx+{u) = Sx_(w) = Sx{oj) and S Y+ {u) = Sy_ M = 
Sy(w)- The phase <j){uj) in the Fig. 5 is chosen to 
achieve the largest possible squeezing, namely, by setting 
dSx{w)/doj = which leads to: 



cos20(w) 



7jg(cjE + n 1 Ki(fc)) 



Q 2 (o;)+7|E 2 (fc) + (n^(fe)7 S ) 2 

2 _ T?2(Z\ , „,2 



(58) 



where O(w) = ur — E 2 (k) + 7^/4, the eg and E(k) are 

related by EqnH] 

Substituting Eqns l51l54l and l55l into Eq. ([56]) leads to 



Sx(u) = 1- 
S y {oj) = 1 + 



N{uo) + "f n nV d (k) 
^nV d {k) 



(59) 



where Af(u>) ■ 
Obviously: 



^{k) + (fiV d (kh % ) 2 



Sx(lo)S y (lo) = AX±AY± = 1 



(60) 



The results show that for a given in-plane momentum k 
and a given individual photon frequency u> with respect 





(a) 



(b) 



FIG. 6: The energy spectrum and the decay rate of the exciton 
versus in-plane momentum k. The direct exciton is always in 
case (a), the quasi-particles in Fig. 4 are not well defined in 
any momentum scales. However, as argued in section IX, 
due to its large lifetime r ex , the indirect exciton is in case 
(b). When |fe| < k* , the quasi-particles in Fig. 4 are not well 
defined . However, when |fc| S> k* , they are well defined 
quasi-particles. Compare with Fig. 4 which is the case with 



to the chemical potential fj,, there always exists a two 
mode squeezing state which can be decomposed into two 
squeezed states along two normal angles: one squeezed 
along the angle 4>(uj) and the other along the angle <f>(u>) + 
tt/2 in the quadrature phase space (X, Y) as shown in 
Fig.5. 

In the following, we discuss two cases \k\ < k* and 
|fc| > k* respectively. For this purpose, we draw the 
exciton energy E(k) and the decay rate 7g/2 in the same 

plot in Fig. 6. When |fe| < k* , the excitons decay very fast 
into photons, so they are not well defined quasi-particles. 
However, when |fc| ^> fc*, the excitons decay into photons 
very slowly, so they are well defined quasi-particles. 

(1) Strong coupling case \k\ < k* : E(k) < 7g/2. 

From Eqnl59( we can see that the maximum squeezing 
happens at a-> m j n = which means at u>k = M : 



Sx(k,Lu = Q) 
cos20(fc, uu = 0) 



i 2 1% iiV d (k) 
N{Q) + nV d (k) 1% 



JV(0) 



(61) 



where W(0) 



Af(uj 



0) 



[E 2 (k) + 7?/4] 2 + (nVd(A;)7j) 2 which is defined below 
Eqnj59l In sharp contrast to the weak coupling case 
E(k) > 7^/2 to be discussed in the following, the 
resonance position u>k = (X is independent of the value 
of fc, this is because the quasiparticle is not even well 
defined in the strong coupling case. The u> dependence 
of Sx{oj) in Eqn l59l is drawn in Fig. 7. The line width of 
the single peak in Fig. 7 is: 



5i(k) 



(62) 



where % = J 4Af(0)[Af(0) + nV d {k) 1% ] 



l\E 2 (k). 



12 



0.2- 
0.0 ■ 



f Ek=0.005pev: 



1.0 
0.8 
_ 0.6 

3, 

W OA 
0.2 
0.0 



-0.02 - 0.01 0.00 0.01 

10 (mev) 
(a) 









/ Ek=0.05(jev 



to (mev) 
(b) 



FIG. 7: The squeezing spectrum at a given in-plane momen- 
tum k when E(k) < 7g/2. The nVd(k) and 7j/2 are fixed at 
50fieV and O.lfieV respectively as determined from experi- 
mental data in section IX. They are also used in the following 
figures 9,10,12-16. There exists only one minimal when the 
photon frequency equals to the chemical potential with the 
width Si (k) given by Eqn l62l Near the resonance, the squeez- 
ing ratio is so close to zero that it can not be distinguished 
in the figure. There are very little difference between the two 
energies E(k) = 0.005^ieV and E(k) = 0.05fieV. This is be- 
cause the quasi-particles are not well defined in this case, the 
squeezing spectrum is mainly determined by nVd(k) and 7j/2. 



(2) Weak coupling case k > k* : E(k) > y%/2. 

As shown in Fig. 6, in this case, |fe| > k* . From EqnloT)! 
we can see that the maximum squeezing happens at the 
resonance frequency w m j n = ±[E 2 (k) — 7§/4] ^ 2 . Recall 
that the Fourier transformation of the Eq. (|4"3"]) gives 
lu = LUk — /i, so the resonance condition becomes u>k = 



H±[E 2 {k)~ 1 l/A] 1 / 2 where 



Sx(k,u; nlin ) = (^i-) 2 = 



hk 2 



E{k)' h 2 k 2 +4MnV d {k) 
cos 2cf)(k, w min ) = 1 (63) 

In this case, (j>(k , w m i n ) = 0, so the two squeezing modes 
are mode + and mode - respectively shown in Fig. 5. 
In sharp contrast to the strong coupling case discussed 
above, the resonance positions depend on the value of 
k, this is because the quasiparticle is well defined in the 
weak coupling case. The squeezing ratio is independent 
of 7£ at the resonances ! Of course, away from the res- 
onances, it will always depends on 7c. From Eqn[S3] we 
can see that increasing the exciton mass, the density and 
the exciton dipole interaction will all benefit the squeez- 
ing. 

The ijj dependence of Sx(w) in Eqn l59l is drawn in 
Fig.8. Whcn£;(fc) > (Qg+^/l + Qi)l%/2, the line width 
of the each peak in Fig.8 is 



S 2 (k) 



(64) 




(a) (mev) 

(b) 



FIG. 8: The squeezing spectrum at a given in-plane momen- 
tum k when E(k) > 7g/2. There exist two minima in the 
spectrum when the photon frequency resonate with the well 
defined quasi-particles. Near the resonance, the squeezing ra- 
tio is so close to zero that it can not be distinguished in the 
figure. When E(k) = 2fieV, the two peaks are still not clearly 
separated. When E(k) — 8/j.eV 2> 7g/2, the quasi-particles 
are well defined which lead to the two well defined resonances 
with width S2(k) in Eqn l65l 



where Qr 



4ftV d (fc) 



It is easy to sec that 82 



JkQk which is equal to the exciton decay rate 7^ multi- 
plied by a prefactor Q^. 

When E(k) < (Qj + + Q k ~)j k '/2, the two peaks are 
too close to be distinguished. 

In short, for a given in-plane momentum, there always 
exists a two mode squeezed state. When E(k) < 7^/2, 
the squeezing spectrum reaches its minimum Eq. [HU 
at u>k = fJ- and the squeezed angle is always non zero 
4>(lo) ^= 0. On the other hand, when E(k) > 7^/2, 
the squeezing spectrum reaches its minimum Eq. [S3] 
at u>k — [E 2 {k) — 7^/4] and the squeezed an- 
gle 0(w m i n ) = 0. In sharp contrast to the widths in 
the ARPS and EDC in Fig. 11, 13 to be discussed in the 
Sec. VII which depend only on 7^, the two widths in 
Eqnl62l and [65] in the squeezing spectra also depend on 
the interaction! The angle dependence of Eqn[5H]in both 
the strong coupling E(k) < 7g/2 and the weak coupling 
E(k) > 7^/2 cases are drawn in the same plot Fig. 9 for 
comparison. From Eq. IrTTl which is valid at E(k) < 7^/2 
( \k\ < k* ) and Eq. [Ml which is valid at E(k) > 7^/2 ( 

I As I > k* ) at the resonance, we can find that the squeez- 
ing ratio and the angle dependence at the resonance on 
the whole in-plane momentum k regime. 

(3) Phase sensitive Homodyne measurement to mea- 
sure the squeezing spectrum and the rotating phase (f>(u>) 

Usual measurements are just intensity measurement 
such as power spectrum experiment and intensity- 
intensity correlation measurement such as HBT 
experimen t 47 ' 48 , so contain no phase information. De- 
tection of squeezed states, on the other hand, requires a 
phase sensitive scheme that measures the variance of a 
quadrature of the photon field. This can be achieved the 
phase sensitive homodyne detection. The experimental 
set-up of this kind of experiment was explained in detail 
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uj (mev) 

FIG. 9: The squeezing angle dependence on the frequency 
when E(k) < 7s /2 and E(k) > j s /2. When E(k) < 7g/2 the 
squeezing angle is always non-zero. Near the resonance, the 
angle is so close to zero that it can not be distinguished in 
the figure. Only when E(k) > 7g/2 and the photon frequency 
resonate with the well defined quasi-particles, the squeezing 
angle is zero. Away from the resonance, the angle becomes 
negative. 



i n 47 i 48 , here we just briefly explain the main points of 
a single mode phase sensitive homodyne detection by a 
schematic Fig. 10. This figure need to be connected with 
the homodyne outputs of the Fig. 16 to detect the two 
modes squeezing spectrum Eqn l59l and squeezing angle 
Eqn l58l It is essentially a phase interference experiment 
between the input light beam A and a local oscillator 
(LO) beam B which is used as a reference beam. Both 
input beam A and LO beam B are incident on a beam 
splitter and are reflected and transmitted, there is a ir/2 
phase shift between the reflected and the transmitted 
beam, so there are two beams C and D coming out 
the splitter, both C and D are linear combination of 
A and B. If one fixed the LO beam B to be a strong 
coherent field with phase <j>, a balanced detection is to 
use a 50/50 beam splitter, the output signal detected by 
the coincidence measurement in the Fig. 10 is taken to 
be the difference between the counting of C photons and 
that of D photons, so it is just the interference between 
the quadrature of the input beam A and the strong LO 
beam B subject to a rotation by angle (j> + tt/2. The 
variance of the output can also be measured which is 
just the squeezing spectrum Eqnl56l By tuning the 
angle 0, both X quadrature and Y quadrature or its 
any linear combination quadrature of the input beam 
can be measured. Then the rotated phase 4> m Eqnl551 
shown in Fig. 5 and drawn in Fig. 9 is just this phase 4> of 
the local oscillator, so they are completely experimental 
measurable quantities in the phase sensitive homodyne 
experiments. 



Input beam 



50/50 
Beam splitter 



Detector 2 



Local Oscillator 



FIG. 10: The balanced homodyne experiment to measure a 
single mode squeezing spectrum. 



VII. ONE PHOTON CORRELATION 
FUNCTION, POWER SPECTRUM AND 
MACROSCOPIC SUPER-RADIANCE 

The one photon correlation function of the output 
field is G±(t) = (a™p(t + T)a°, u l{t)) and the angle 

\ ±k ±k I in 

resolved power spectrum (ARPS) of the output field 
is S±(k,uj) = /+~ dre- iUT G±(T). The normalized 
first order correlation function is defined by g±(r) — 
G±(t)/G±(0). In the following, we will evaluate these 
quantities respectively. 

(1 ) The angle resolved power spectrum (ARPS) 
By inserting Eqn l51[ one obtains the photon number 
spectrum S±(k,uj) — Si(k,uj) and 



= -(S x (w) + Sy(w)-2) 
= 7|» 2 ^ 2 (fc) 



-7p 2 (fc) 



(65) 



7^/4. From Eqnl6"01 one can 



where Q(ui) = uj 2 — E 2 (k) 
see that Sx (uj) + Sy(uj) > 1, so the more squeezing, the 
stronger the power spectrum. If there is no squeezing 
Sx(uj) — S'y(w) = 1, then there is no power emitted, 
this is just the input vacuum state in the Fig. 4 which is 
a coherent state itself. The angle resolved power spectra 
(ARPS) with different E(k) and 7 £ are shown in Fig. 11. 
The total ARPS is the sum of the condensate and the 
quasi-particles: S(k,uj) = So(k,uj) + Si(k, uj). 

In the strong coupling case k < k*,E(k) < 
7g/2, Si(k, us) reaches the maximum •j 2 .n 2 V^(k)/[ r y 2 /4: + 



E 2 {k)} 2 
Si(k,uj) 



0, E(k) = u\k\ -> 0, then 
so the curve has a half width 



at ujk = A*- As k 

7|» 2 V|(g) 
(w 2 + 7 H/4) 2 ' 

~ ?l7o ~ 10~ 4 mey . This is expected, because the quasi- 
particles are not well defined with the decay rate 70 much 
larger than its energy E{k). 

In the weak coupling case E(k) > 7g/2, at the two res- 
onance frequencies uj k = fi± [E 2 (k) — 7S/4] 1 / 2 , Si(k,uj) 

reaches the maximum n 2 V^ (k) / E 2 (k) which only de- 
pends on the exciton density, the dipole-dipole interac- 
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E*=0.05 Mev 
Ek=0.5 Mev 
Ek=2 Mev 




S(k )/L z 



FIG. 11: The angle resolved power spectrum (ARPS) of the 
emitted photon with in-plane momentum k. When E(k) < 
hfj:/2, there is only one peak in the power spectrum with 
the width 7g/2. When E(k) > ft7g/2, there are two peaks 
in the power spectrum at the two resonance photon frequen- 
cies uik — M i E(k) also with width 7g/2. The Si(k,u>) at 
E(k) = 0.5/zeV and E(k) = 2^eV are multiplied by 50 and 
500 in order to be seen in the figure. Compared to the squeez- 
ing spectrum in Fig. 9, one can see the ARPS can distinguish 
the two quasi-particle peaks clearly even at E(k) — 0.5fieV . 
While the squeezing spectrum in Fig. 9 can not distinguish the 
two quasi-particle peaks even at E(k) = 2/j.eV, so the ARPS 
is a much sensitive probe of the quasi-particle spectrum than 
the squeezing spectrum. 



tion and the quasi-particle spectrum, but independent of 
7g ! It can be shown that when E(k) ^ 7^/2, the width 
of the two peaks at the two resonance frequencies is ~ 7g, 
this is expected, because the quasi-particle is well defined 
with energy [E 2 (k) — 7S/4] 1 / 2 and the half-width 7^. In 
sharp contrast to the two widths in Eqn l62l and l65l in the 
squeezing spectra which depend on both the interaction 
and 7g, the widths in the ARPS and EDC in Fig.11,13 
depend only on jt- 

(2) Momentum Distribution Curve (MDC) 
The power spectrum at a given in-plane momen- 
tum k is Si(fc) = J2k z Si(k,u) = J duo k D^(uo k )Si{k,ijj) 
which is nothing but the Momentum Distribution Curve 
(MDC^i: 



Si(k) 



oc L, 



(66) 



As to be explained in section IX and shown in the Fig. 6, 
the crossing point is at k* ~ 10~ 2 em _1 . From Eqnl3UI 
one can see the condensate contribution at k = is 
N p h/L z = Njo/v C j oc N, while the contribution from 

the quasi-particle is Si(k — > 0)/L z — 2 " t) V ^J°^ oc n 2 /7o- 
So the MDC is a Lorentian with the half width at k* 
as shown in the Fig. 12. So the k* has a clear physical 
meaning as the half width of the MDC and is an experi- 
mentally measurable quantity. 

In fact, we can also calculate the one photon correla- 




FIG. 12: The zero temperature Momentum Distribution 
Curve (MDC) has a bi-model structure consisting of a sharp 
peak S(k — 0)/L z = N^yo/vg due to the condensate at 
k — superposing on a Lorentzian peak with a half width 
k* ~ 10 2 cm -1 due to quasi-particle excitations at k 7^ 0. The 
coherence £ ~ 1/fc* has been measured in 2 ^. 



tion function: 



G(r) 



d 2 k 



(2tt) 2 fc 2 + fc* 2 



(67) 



where we can identify the coherence length £ ~ 1/fc* ~ 
40/xm. This coherence length has been measured inSS, 
and will be discussed in detail in section IX. 

(3) Energy Distribution Curve (EDC) 

The power spectrum at a given energy u> is Si(uj) = 
^2j;Si(k,u)) which is nothing but the Energy Distribu- 
tion Curve (EDCp-: 



Si(w) = Nnx 



d 2 k 



(2?r) 2 + 1 2 E 2 (k) 



NhV$(k = 0) M 



f{^)cxNn 

Tfc=0 



(68) 



where f(x) = — arctg 1 ^ 4 X X } where — 7r/2 < 

arctgy < tt/2. Because /(0) = 4, then Si(w = 0)/N = 



nV d 2 (fc=0) 



The S\(w) is shown in Fig. 13a. 



The Si(u) is the sum over all the angle resolved power 
spectrum curves in Fig. 11 at the fixed energy lj = Lo k — [i. 

As fc - 0, E(k) = u\k\ -> 0, then Si(k,u) -> gr+^r , 

so the curve has a half width ~ S.70 ~ lO^^meV. As 
|fc| increases to fc* ~ 10 _2 cro -1 , the curve starts to 
split into two peaks as shown in the Fig.ll, then when 
\k\ — > |fc OT aa;|/2, the two peaks stand for the two well de- 
fined quasi-particle excitations, the splitting of the two 
peaks reaches ~ u\k max \ ~ O.lmeF. So the EDC curve 
will simply smear out all the fine structures of the an- 
gle resolved power spectrum in Fig.ll and finally end up 
with an envelop curve with a half width w e ~ 7^ as shown 
in the Fig. 13. From Eqnl37l one can see the condensate 
contribution at fc = is n Wk /N — 7o<5(wfc_ — fi), while 
the contribution from the quasi-particle S%(uj = 0)/N = 
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S(f> )/N 




Quasi-particle^^ 



W e W„ 
H=1.54eV 

(b) 




FIG. 14: The one photon correlation function 



FIG. 13: (a) The Energy Distribution Curve ( EDC ) from the 
quasi-particle. (b) The zero temperature EDC has a bi-model 
structure consisting of a sharp 8 function peak S(uj)/N — 
7o<5(w — fj) due to the condensate at k — superposing on 
a Lorentzian peak with the half envelop width w e ~ 7g ~ 
0.1 jj,eV due to quasi-particle excitations at k ^ 0. 
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i — as shown in the Fig. 13a. As to be explained in 
section IX, the EDC is what experiments ini£ measured 
at various P ex , V g and T. 

(4) The total radiation rate from the quasi-particles 
The total number of photons emitted from the quasi- 
particles is: 



Sr = J2J2 Sl $,u)~NxL z 



(69) 



which is proportional to the total normalization volume 
of the system as it is expected. This can also be used 
as a self-consistency check on our results achieved on the 
quasi-particle part at k ^ 0. As shown in the section IV, 
this self-consistency check played very important roles on 
the condensate part k = 0. In fact, if taking Eqns |6"51 and 
l6"6"lat face value, then the total number of photons Eqnl6"9l 
diverge, but this should not cause any concern, because 
both Eqnsl6"51 and 1661 only hold at small momentum k < 
k m ax and to < uk max respectively. 

The radiation rate along a given direction (k, k z ) is: 



P± (^j 



L 2 L Z 



oj k v„ x L 



(70) 



which vanishes in the limit L z — > oo as expected. 
Then the radiation rate at a given in-plane momentum 
P[ a (k) = VgfiSi(k)/ L z and the radiation rate at a given 
energy Pf a (uj) — v g u,S\(w)/L z where we have used the 
fact that S\{k,ui) is a even function of w = lo^ — /■*• 
The total radiation rate from the quasi-particles is: 

p r = Prd & = VV £ * cx v g u.Nn (71) 



where again k z is sum over both the upper and the lower 
space in the Fig. lb. P[ a is also cx fiN just as the ra- 
diation rate P™ from the condensate. Because of the 



weak dipole-dipole interaction, the quantum depletion is 
small, so P[ a is still much smaller than Pg° in Eqnl38l 
Furthermore P[ a is spread over all the possible angles, 
while P™ is focused along a single direction. 

(5) One photon correlation functions 

By using the Fourier transformation to Si(k,u), one 
can get G±(k, r) = G\(k, r) 



Gi(fc,r) 



4P(fc) ^ E(k) 



+ h.c], (72) 



The normalized first order correlation function 
g±(k,r) =pi(fc,r), where 



g\{k,r)=e — T [cos(E(k)T) + 



Ik 



2E(k) 



sin(P(fc)T)] (73) 



It turns out that the first order correlation function is 
independent of the relation between E{k) and 7^/2 and is 
shown in Fig. 14. The G\ was measured in the EHBL in 25 
at T = 1.6K and in exciton polariton ir>22 at T = 4K. 
The effects of finite temperature and trap potential must 
be considered before comparing our theoretical results 
with the experimental data. 

(6) The quasi-particle spectrum in a non-equilibrium 
stationary exciton superhuid 

It is important to compare the excitation spectrum 
in Fig. 4, Fig. 6a and Fig. 6b. Fig. 4 is the well know 
quasi-particle excitations in an equilibrium superfluid. 
They are well defined quasi-particles with infinite life- 
time. However, the quasi-particles in Fig. 6a are not well 
defined in any length scales, because the decay rate is 
always much larger than the energy. Fig. 6(b) is between 
the two extreme cases. When k < k*, the quasi-particle 
is not well defined, the ARPS is centered around LOk = \i 
with the width 7^. The MDC has large values at k < k*. 
The EDC has large values at lu < 7^. When k > k* , the 
quasi-particles is well defined, the ARPS has two well de- 
fined quasi-particles peaks at ujk = /U± [E 2 (k) — 7I/4] 1 / 2 
with the width 7^. The MDC has very small values at 
k > k* . The EDC has very small values at u> > 7&. So 
in the long wavelength r > £ ~ 1/k* ( or small momen- 
tum k < k* ) limit and long time r > r ex ~ l/7fc (or low 
energy limit lu < 7^ ) limit, there is not a well defined su- 
perfluid which is consistent with the results achived in 52 . 
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However, in the distance sacle r < £ ~ 1/fc* (or momen- 
tum k > k* ) limit and the time scale r < r ex ~ 1/7^ (or 
energy scale w > 7^ ) , there is still defined superfluid and 
associated quasi-particle excitations. This is the main 
difference and analogy between the equilibrium super- 
fluid in Fig. 4 and the non-equilibrium steady state super- 
fluid in Fig. 6b. Very recently, the elementary excitation 
spectrum of exciton-polariton inside a micro-cavity was 
measured 3 -^ and was found to be very similar to that in 
a helium 4 superfluid shown in Fig. 4 except in a small 
regime near k = 0. We believe this observation is pre- 
cisely due to the excitation spectrum in a non-equilibrium 
stationary superfluid shown in Fig. 6b. 

(7) The Superradiance from the quasi-particles 

Note that the angle resolved power spectrum, the MDC 
and EDC in Eqns l6"5"ll6f)l |6"51 are all proportional to N 2 
instead of N. It is the characteristic of super-radiance 
in a macroscopic system. This should not be too sur- 
prising, because the excitonic superfluid is a macroscopic 
quantum coherence phenomena, so it is natural to lead 
to macroscopic superradiance. As k — > in the Fig. 6, 
Sl (k) ~ ^M. ~ n 2 / 1% where 7/c appears in the de- 
nominator in the strong coupling case, so the macro- 
scopic superradiance can only be achieved by the non- 
perturbative calculations presented in this section, but 
can not be derived by any finite order perturbative cal- 
culations presented in the appendix D. 

In conventional quantum optics, N two level atoms 
interacting with a single ( or multi-) photon mode(s) in- 
side a cavity. If the N static atoms are confined into a 
small volume V inside the cavity which is much smaller 
than the wavelength of the photon mode, then the in- 
teraction between all the N atoms and the photon mode 
can be taken as the same constant A, then when half of 
the atoms are in the excited level, the radiation intensity 
from the N atoms is proportional to N 2 instead of just 
N during the time interval ~ 1/-ZV, so the total power 
emitted during this time period is N 2 x 1/N ~ N as 
required by the energy conservation. This is due to the 
cooperative effects of the N atoms which is due to the 
fact that the N atoms, being interacting with the same 
photon field, so can not be treated as independent N 
atoms. This is called non-equilibrium superradiance first 
studied by Dicke^. Generalizing the Dicke model to a 
stationary state inside a cavity was studied i n 69 i 70 . It 
was found that there is a second order phase transition 
driving the coupling constant A: when A < A c , the sys- 
tem is in a normal phase, when A > A c , the system is 
in a superradiative phase. It was also pointed out inZS 
that it is very unrealistic to realize Dicke model in the 
thermodynamic limit N, V — > 00, but keep N/V to be 
finite, because it is essentially impossible to make V still 
smaller than the wavelength of the photon in the ther- 
modynamic limit. So the superradiance is essentially an 
effect for finite number of static atoms confined into a 
small volume. 

The superradiance from the ESF has completely differ- 
ent mechanism: (1) the size of the sample is much larger 



than the wavelength of the photon field (2) The photons 
field is a continuum of photons with different k z at a 
given k in Eqn ll5l so acting as a reservoir to the excitons 
with in-plane momentum k (3) all the excitons are al- 
ways in motion, in fact, moving in a coherent fashion. So 
all these conditions violate the conditions to achieve the 
superradiance in conventional quantum optics. So the 
collective radiation from the ESF is due to the macro- 
scopic coherence of the exciton superfluid itself which is, 
in turn, due to the dipole-dipole repulsion. 



VIII. TWO PHOTON CORRELATION 
FUNCTIONS AND PHOTON STATISTICS 

The quantum statistic properties of emitted photons 
can be extracted from two photon correlation functions. 
The normalized second order correlation functions of the 
output field for the two modes at k and — k are 



and 



92 ( T ) 



out\ 

a- 

k 



{t)a°£ t \t + T)aft{t + T)a°f{t)) i 



|Gi(0)r 



(74) 



(75) 

The second order correlation function (r) deter- 

mines the probability of detecting n_£ photons with mo- 
mentum — k at time t and detecting photons with 

momentum k at time t + r. Just like the one photon 
correlation function in Eqnl73[ it turns out that the sec- 
ond correlation functions are also independent of the re- 
lation between E(k) and 7^/2 and are shown in Fig. 15. 

The Wick theorem^ gives g£' '(r) = 1 + |3i(t)| 2 where 
gi(r) is given by Eqn(73] Similarly, it can be shown that 

g { f k \r) = 1 + l/i (r)| 2 where 

E 2 (k) + 4 -1£.t 
/i(t) = e 



x \ul 



nV d (k) 

e -iE(%)i 



iE{k)n 



+ v z r - 



L ' T - li: ' k E(k)+i^- 



(76) 



'E{k)-if 

Then the normalized second order correlation functions 



gf\r) = 1 + e-^[cos(E(k)r) + - * sm(E(k)r)} 2 

(77) 

and 



[cos(E(k)T) 



7fe 



2E(k) 



sin(E(k)r)] 2 l78) 



17 



M 
2.5 
2D 

(±10 

g 2 (x) 15 

Lfl 
0.5 
(111 



homodyne 



homodyne 




x (ns) 



FIG. 15: The two photon correlation functions between k and 
— k against the delay time r. The envelope of second order 
correlation function decreases as time evolves indicates that 
the emitted photons is photon bunching and the photo-count 
statistics is super-Poissonian. 



It is easy to see that the envelope decaying function 
is given by the exciton decay rate 7c, while the oscil- 
lation within the envelope function is given by the Bo- 
goliubov quasi-particle energy E(k). Subtracting EqnlTSl 
from Eqnj77]lead to: 



(±k)t \ (fc), , E 2 (k) 
92 V ) 92 ( r ) = e 
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4 



(79) 



So we only draw (t) in the Fig. 18. 

When r = the two photon correlation function are 

92°\ty ~ 2, so just the mode k alone behaves like a 
chaotic light. This is expected because the entanglement 
is only between — k and k. In fact, 



# fc) (0) 



E 2 (k) + ? 



> .9^(0) = 2. 



(80) 



So it violates the classical Cauchy-Schwarz inequality 
which is completely due to the quantum nature of the 
two mode squeezing between k and —k. 

The normalized two photon correlation function 

92 ( r ) i s shown in Fig. 15. From Fig. 15, we can 
find that the two photon correlation functions decrease 
as time interval r increases which suggests quantum na- 
ture of the emitted photons is photon bunching and the 
photo-count statistics is super-Poissonian. See Fig. 17 for 
its experimental measurement. 



IX. DISCUSSIONS ON AVAILABLE 
EXPERIMENTAL DATA AND POSSIBLE 
FUTURE EXPERIMENTS 

We will determine how all the important parameters in 
our theory can be precisely measured by the experiments. 
Then we will compare our results on the EDC with the 




Coherent light 

FIG. 16: The experimental set-up for the homodyne detec- 
tion of the two modes squeezed photon in the EHBL system. 
The coherent photons are emitted along the z-direction re- 
flecting the nature of the condensation. The squeezed pho- 
tons are emitted along all the tilted directions reflecting the 
fluctuation above the condensation. There is a macroscopic 
superradiance when the angle 9 is small. The homodyne com- 
ing out of or A°y- l _ is connected to Fig. 10 for the single 
mode squeezing spectrum measurement. 



available experimental data^ and then discuss possible 
future experimental set-ups to detect all the other theo- 
retical predictions achieved in the previous sections. 

1. Comments on current experimental data 

Ini£, the spatially and spectrally resolved photolumi- 
nescence intensity has a sharp peak at the emitted photon 
energy E — \.bAheV with a width ~ ImeV at the low- 
est temperature ~ IK ~ O.lmeV. The radiation power 
P ex varies from YlpW to 1.6mW. The gate voltage is 
set at V g = 1.2V. As shown in Eqn[T|)l one can iden- 
tify E = (i = E g - E b + nV d (0) = 1.545eV. The energy 
conservation at k z = gives the maximum in-plane mo- 
mentum k max ~ fi/vg ~ 4.3 x 10 4 cm _1 where we used 
the speed of the light v g ~ 8.7 x 10 9 cm / s in GaAs^. Then 



the maximum exciton energy E„ 



Q.lhmeV 



where we used the spin wave velocity u~5x 10 5 cra/s. 
Taking the mass of the exciton M ~ Q.22m^, then from 
the expression of the spin wave velocity u — yJnVd{fy/M 
calculated in EqnJ^H we find nVd(0) — 0.05meV which 
is the value we used in all the previous Figs. 7-15. Tak- 
ing the exciton density n ~ 10 10 cto~ 2 , we can see that 
K = 2ir/a ~ 10 6 cm _1 where a is the average spacing 
between excitons. The average lifetime of the indirect 



excitons in the EHBL 4 '- is r, 



40ns, then we can esti- 
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mate the exciton decay rate 7fc ~ l/40ns ~ 10~ 4 meV = 
O.lfxev. At the boundary of the two regimes in the Fig. 6 
where E(k*) — uk* — 7^* /2 = 10~ 4 meV, we can extract 
k* = 2.4 x 10 2 cto _1 . So there are three widely separated 



of the photo-detector in the Fig. lb, one should see 



momentum scales, k* ~ 10 2 cm 1 -C fc„ 



10 4 cm- 1 < 



if ~ 10 6 cm 1 ~ fe r which is the roton minimum in Fig. 4. 
So all the important parameters such as the chemical po- 
tential u, the quasi-particle energy E(k) — u\k\, the ex- 
citon decay rate 7^ and the exciton dipole-dipole inter- 
action strength nVd(0) in our theory can all be extracted 
from experimental data. 

The typical trap size ( or exciton cloud size ) L ~ 
30/ito. The number of excitons is N = nL 2 ~ 10 5 which 
is comparable to the number of cold atoms inside a trap 
in most cold atom experiments. The central peak due 
to the condensate in the MDC Fig. 12 is broadened to 
ko ~ 1/L ~ 10 3 cto _1 which is already larger than the 



half width due to the quasi-particle k* ~ 10 



So 



it is impossible to distinguish the bi-model structure in 
the MDC at such a small exciton size. The coherence 
length was measured in^ 5 .. From Eqn l67l we find the 
coherence £ ~ 1/k* ~ 40/xm which is slightly larger than 
the exciton cloud size L ~ 30/im. 

For an inhomogeneous condensate (b(k)) — ip(k) inside 



a harmonic trap V(r) 



we assume local density 



approximation (LDA) is valid, then Eqn l26l should be 
replaced by: 

idta(k, k z ) — (uj(k,k z )—fi+V(r)—iK)a(k 1 k z )—ig*(k,k z )'ip(k) 

(81) 

where we have still dropped the zero mode fluctuation 

of the inhomogeneous condensate. Because g*(k,k z ) ~ 
1/2 

L z 1 -► and k = v g / L z — > such that a stationary 
state is reached where the energy u)(k,k z ) is pinned at the 
local chemical potential /i(r) — \i — V(r), so the central 
peak due to the condensate in the EDC in the Fig. 13a is 
broadened simply due to the change of the local chemical 
potential from the center to the edge of the trap Aui — 
O.lmeV where we used the experimental valued 
10~ 12 eVnm~ 2 . This geometrical broadening due to 



\ul? 



the trap is already much larger than the half width due 
to the quasi-particle w e ~ 7^ ~ O.lfieV in Fig. 13b, so 
it is impossible to distiguish the bi-modal structure in 
the EDC in the Fig. 13b. In order to understand if the 
observed EDC peak with the width ~ lmeV at the lowest 
temperature ~ 1.7 K in the experiments in 4 ^ is indeed 
due to the exciton condensate, one has to study how the 
bi-model structures shown in Fig. 13b will change inside 
a harmonic trap at a finite temperature ~ 1.7 K and the 
effects of both dark and bright excitons. This will be 
discussed in a separate publication 53 . 

2. Possible future angle resolved experiments 
As shown in section VII, the MDC, especially the EDC 
smear out all the fine structures of the angle resolved 
power spectrum in the Fig. 11. So in order to test the 
theoretical results precisely, it is necessary to perform the 
ARPS measurement. As one rotates the angle tan0& = 



the following interesting behaviors. When \k\ < \k*\, the 
ARPS and the squeezing spectrum only have one peak 
at the resonance frequency u>k = u ~ 1.545eV as shown 
in Fig. 11, Fig. 7 and Fig. 18b. The line width of the peak 
is uniquely determined by the quasi-particle excitation 

k and the decay rate 7?. The one peak will start to 

split into the two peaks at u k* =7g,/2 corresponding 
to the angle sin#* = k*/k rnax ~ 10 -2 , so 9* ~ 10~ 2 . 
When k>k*( namely, 9 > 9* ), both the ARPS and the 
squeezing spectrum have two peaks at the two resonance 

frequencies uu k = v g (k 2 z + fc 2 ) 1 / 2 = 1.545eV± [u 2 k 

72/4] 1 / 2 as shown in Fig. 11, Fig.8 and and Fig.l8b. The 
position and the width of the two peaks are uniquely 
determined by the quasi-particle excitation uk and the 
decay rate 7^. So all the characters of the condensate 
and the fluctuations above it are reflected in the angle 
resolved measurements of squeezing spectrum and the 
power spectrum. 

3. Possible future two modes phase sensitive homodyne 
experiments 

As shown in the appendix C, due to the relation 
Eqn lC5[ the experimental set-up in Fig. 16 combined 
with the single mode phase sensitive homedyne set-up 
in Fig. 10 can measure the squeezing spectrum in Eqnl59l 
and the squeezing angle in Eqnl58l 

4. Possible future two modes HanburyBrown-Twiss 
type of experiments 

The single mode two photon correlation functions 

92 ( T ) m Eqnj74] can be measured by the usual HBT 
set-up 47 ' 48 . It is not very interesting anyway. From the 
relation in Eqn ICli we can see the experimental set-up 
in Fig. 17 can measure the most interesting two modes 

two photon correlation functions g^^M m Fan. 1781 and 
Fig.15. 



X. CONCLUSIONS 

In conventional quantum optics, coherent light was 
produced by stimulated radiations from a pump which 
leads to particle number inversion and amplified by opti- 
cal resonator, here in EHBL, the coherent state along the 
normal direction is due to a complete different and new 
mechanism: spontaneous symmetry breaking. It has the 
following remarkable properties: (1) highly directional: 
along the normal direction (1) highly monochromatic: 
pinned at a single energy given by the chemical potential 
/i (3) high power: proportional to the total number of 
excitons. These remarkable properties could be useful to 
build highly powerful opto-electronic device. 

In conventional non-linear quantum optics, the gen- 
eration of squeezed lights also requires an action of a 
strong classical pump and a large non-linear suscepti- 
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Coherent slate 



FIG. 17: HanburyBrown-Twiss ( HBT ) type of experiment 
to measure the two modes two photon correlation functions 
^ ±fe) (r) in Eqn[78]and Fig. 15. 

bility ■ The first observation of squeezed lights was 
achieved in non-degenerate four-wave mixing in atomic 
sodium in 1985 5 -. Here in EHBL, the generation of the 
two mode squeezed photon is due to a complete different 
and new mechanism: the anomalous Green function of 
Bogoliubov quasiparticle which is non-zero only in the 
excitonic superfluid state. The results achieved in this 
paper are robust against any microscopic details. The 
applications of the squeezed state include (1) the very 
high precision measurement by using the quadrature with 
reduced quantum fluctuations such as the X quadrature 
in the Fig. 2 and 3 where the squeeze factor reaches very 
close to at the resonances (2) the non-local quantum 
entanglement between the two twin photons at k and — k 
can be useful for many quantum information processes. 
(3) detection of possible gravitational waves^i. 

In conventional quantum optics, N non-interacting 
atoms confined in a small volume interact with photon 
modes, a super-radiance was proposed by Dicke. Very re- 
cently, there are preliminary experimental evidence that 
excitons in assembly of quantum dots may emit Dicke's 
superradianceZ^. In conventional laser, due to random 
motions of atoms confined in a volume larger than the 
wavelength of the laser beam, the laser power is only pro- 
portional to N when the pump is above the threshold. So 
far, no super-radiant laser has been achieved. Here the 
ESF phase is a macroscopic quantum phenomena in a 
macroscopic sample ( namely in the thermodynamic limit 
), so the super-radiance emitted from this system has 
completely different mechanism than the Dicke model. 

In conventional quantum phase transitions, all the 
quantum phases and phase transitions are stable phases 
at equilibrium. For example, all the possible interesting 
quantum phases in BLQH mentioned in the introduction 
are stable equilibrium phases. However, the quantum 
phases in EHBL in Fig. 3 are just meat-stable phases, they 
will eventually decay through emitting photons, photons 



are very natural internal probe of the quantum phases 
and quantum phase transitions. The characteristics of 
the photons such as power spectrum, squeezing spectrum 
and photon correlations are completely determined by 
the nature of the quantum phases such as the ground 
state and elementary excitations. The excitation spec- 
trum in a non-equilibrium superfluid shown in Fig. 6b is 
also different from that in a conventional equilibrium su- 
perfluid shown in Fig. 6a. This difference completely and 
precisely explained the recent experimental observation 
of excitation spectrum of exciton-polariton in a planar 
microcavity^l. 

In conventional condensed matter experiments, ground 
states are completely stable, so in order to externally 
probe the quasi-particle excitations of a quantum phases, 
the temperature to perform the experiments has to be 
sufficiently high so there are enough quasi-particle exci- 
tations excited above the ground state. However, the 
quantum phases in the Fig. 3 are just meta-stable, so 
the internal probe of emitted photons from the quantum 
phase can also reflect the energy spectrum of the quasi- 
particle even at T = 0. This is the salient feature of the 
internal probe in meta-stable systems different from the 
external probes in conventional stable condensed mat- 
ter experiments. Exploring the connection between the 
quantum phase transitions and quantum optics in meta- 
stable non-equilibrium systems is still a developing field 
and very exciting and rich. 

It is very instructive to compare the possible exci- 
ton BEC in EHBL with the well established BEC of 
ultra-cold atoms^. The quantum degeneracy temper- 
ature of two dimensional excitons can be estimated to be 
T% x ~ 3K for exciton density n ~ 10 10 cto -2 and effective 
exciton mass m = 0.22m e , so it can be reached easily by 
He refrigerator. While due to the heavy mass of atoms 
and the dilute density, T% tom ~ [iK. Both the exciton 
BEC and the atomic BEC belong to the weakly inter- 
acting BEC class, so Bogoliubov theory apply to both 
cases. The BEC to BCS crossover and quantum phases 
in Fig. 3 tuned by r s or imbalance is also quite similar to 
those of two species of ultra-cold neutral fermions tuned 
by Feshbach resonance or imbalance^. In fact, atomic 
BEC is also a meta-stable ground state, it will eventu- 
ally evaporate away. The conventional way to detect the 
atomic BEC is through the time of flight experiments 
which will destruct the BEC. The smoking gun experi- 
ment to prove the realization of BEC is through the ob- 
servation of vortices when the BEC is under rotation. In 
this paper, we explicitly showed that although it is hard 
to rotate the metastable excitonic BEC to generate the 
superfluid vortices, the non-equilibrium stationary Bo- 
goliubov quasi-particle spectrum of the exciton BEC in 
teh Fig. 6b can be directly extracted from various char- 
acteristics of the emitted photons even at T = 0. A 
superradiant behavior was found in the off-resonant light 
scattering experiment from the BEC condensate^ and 
studied theoretically ir*Z2. In the off-resonant scatter- 
ing, the excited atomic state is adiabatically eliminated, 
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so important phase information containing the excita- 
tion spectrum is lost in the adiabatic elimination proce- 
dure. In the exciton BEC phase studied in this paper, 
the photons are always internally resonant with the exci- 
tons, so it becomes a complete natural internal probe of 
the ground state and excitations of all the possible exci- 
ton quantum phases. However, so far, there is very lit- 
tle experimental ways to measure the Bogoliubov quasi- 
particle of atomic BECS. So detecting quantum phases 
in ultracold atoms remains an outstanding problem. Re- 
cently, one of the authors and his collaborators devel- 
oped a theory to detect the nature of quantum phases of 
ultra-cold atoms loaded on optical lattices by using cavity 
enhanced off-resonant light scattering™ Excitons carry 
a electric dipole moment. While in the context of cold 
atoms, very exciting perspectives have been opened by re- 
cent experiments on cooling and trapping of 52 Cv^- and 
polar molecules 5 ^.. Being electrically or magnetically po- 
larized, the 52 CV atoms or polar molecules interact with 
each other via long-rang anisotropic dipole-dipole inter- 
actions. The superfluid and solid phases of these polar 
molecules are under extensive experimental searc h 58 ' 60 . 
The excitation spectrum very similar to that in Fig. 3 
including the roton minimum has also been proposed in 
trapped pancake dipolar Bose-Einstein condensates^ - in 
atomic experiments. The presence, position and depth of 
the roton minimum are tunable by varying the density, 
confining potential. So the insights and results achieved 
in this paper may also shed some lights on cold atoms 
and molecules. 

Before getting to the final summary of the results 
achieved in this paper, we briefly mention some previ- 
ous and very recent work on semi-conductor electron- 
hole bilayers. BCS pairing of excitons and BEC to BCS 
crossover were discussed i n 62 i 63 . These mean field calcu- 
lations can only describe the quasi-particle part, but con- 
tains no information on the collective modes. However, 
on the BEC side in the Fig. 3 discussed in this paper, it 
is the collective mode shown in Fig. 6b which contributes 
to the two mode squeezing state. Quantum Monte-carlo 
using trial- wavefunctions on the global phase diagram in 
Fig. 3 except the possible excitonic supersolid phased 5 - 
were studied i n 64 ' 65 . The transport properties were stud- 
ied in^ 6 ., but we disagree with the claim made that 
an in-plane magnetic field can induce a counterflow su- 
percurrent. The photoluminescence from the excitons 
inside a harmonic trap was studied by using first order 
Fermi Golden rules in££ by treating the exciton BEC as 
a single two level atom. This Golden rule treatment was 
not able to capture any of the physical phenomena ex- 
plored in this paper which treat the exciton BEC as a 
quantum phase with its symmetry breaking ground state 
and Bogoliubov excitation spectrum. In the appendix D, 
we will use a first and second order Golden rule calcula- 
tions which treat the exciton BEC as a quantum phase 
with a ground state and infinitely many excited states 
with different number of Bogoliubov quasi-particles. This 
kind of Golden calculation based on many particle ground 



and excited states can see some signatures of the results 
on coherent state and squeezed state achieved by non- 
perturbative calculations in the main context. The two 
phonon squeezing generated by lattice an-harmonic ef- 
fects and second-order Raman scattering was discussed 
ir*Z£. As shown in section VII, the squeezing is closely 
related to the macroscopic superradiance, so it would be 
interesting to see if one can achieve a macroscopic super- 
radiance of phonons. 

Several kinds of measurements such as angle resolved 
power spectrum, two mode squeezing spectrum, one pho- 
ton and two photon correlations functions can completely 
detect the characters of the emitted photons which, in 
turn, are very natural internal probes of the ground state 
and excitations of exciton of the quantum phases in the 
EHBL system. We established the direct relation be- 
tween the nature of quantum phases and the nature of 
the emitted photons from the quantum phases: the BEC 
of excitons lead to the coherent state ( or a spontaneous 
laser ) of emitted photons. The anomalous Green func- 
tion of the Bogoliubov quasi-particles lead to the two 
mode squeezing of emitted photons along all the titled 
directions. The whole system behave coherently to emit 
superradiance even in the thermodynamic limit. In fact, 
the ESF phase of the excitons play a similar role as a two 
mode squeezing operator which squeezes the input vac- 
uum state into a two mode squeezed state in a suitable 
rotated frame. We also evaluated the energy distribution 
curve (EDC) and momentum distribution curve (MDC) 
and compared our results with the available experimen- 
tal data on MDC and EDC. We determined how all the 
important parameters such as the chemical potential /i, 
the quasi-particle energy E(k) = u\k\, the exciton decay 
rate 7g and the exciton dipole-dipole interaction strength 
nVd(0) in our theory can all be extracted from experimen- 
tal data. The photons in all directions show bunching and 
super-Poissonian. Possible future angle resolved power 
spectrum experiments, two modes phase sensitive homo- 
dyne experiment to measure the squeezing spectrum, two 
modes HanburyBrown-Twiss experiment to measure the 
two photon correlation functions can be used to directly 
test these predictions. In fact, if the exciton BEC was 
indeed achieved in the past experiments^' 5 -' 6 - remains un- 
clear. In future publications, we will study the effects of 
traps, spins of the excitons and finite temperature. The 
theoretical results achieved and possible future experi- 
mental set-up proposed in this paper should shed consid- 
erable lights to test if the exciton BEC can be observed 
in the EHBL without any ambiguities. 
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APPENDIX A: THE PHYSICAL PICTURES IN 
TERMS OF RADIATION ZONES 

For simplicity, we first explain the physical picture at 
the zeroth order where we neglect the decay of excitons 
into photons in the Fig. 18a. Note that at zero tempera- 
ture T = 0, there is no quasi-particle (3% excitation in 
Eqnl2T)l namely, the input state in the Fig. 4. is the 
ground state \BEC > ! However, what couples to the 




(a) 



(b) 



FIG. 18: The radiative zone (RZ) in BEC ESF phase (a) 
The exciton lifetime T ex — * oo. The RZ is the regime of E(k) 
between the four intersection points, (b) The exciton lifetime 
r ex is large, but finite. The solid gray parabola lines denote 

H ± [u 2 Ifcl 2 - 7g/4] 1/2 when E(k) > 7g/2, the blue dashed 



photons is br in EqnUS] which is a linear combination of line denote w k = v g J \k\ + ki at k z = 0. When E(k) < 



/?£ and pi. Equivalently, the phase 9 fluctuation used in^ 
is a relativistic excitation which has both positive parti- 
cle branch and negative energy anti-particle branch. The 
vacuum energy of the photon E^ h is sitting below that 
of the exciton Eq x by the chemical potential /i, namely, 
Eq x — Eq = /i. There exists the resonance condition not 
only between the photon and the positive energy branch 



of the particle: hu k = 1.545eV- 



, but also between 



the photon and the negative energy branch of the anti- 



particle: fiiij-k, = 1.545eV- 



as shown in Fig. 18a. 



When the effects of the excitons decaying into pho- 
tons are taken into account self-consistently, the above 
picture need to be modified especially at small value of 
k. When |fe| < k* , the two quasi-particle branches in 
Fig. 18a are not well defined anymore, the power spec- 
trum and the squeezing spectrum only have one peak at 
resonance frequency u>j, = fj, = 1.545eV. The one peak 

starts to split into two peaks at u k* =7£» /2 where the 
quasi-particle excitations are still not well defined. How- 
ever, when |fc| 3> |fc*|, the two quasi-particle branches 
in the Fig. 18a remain well defined, so both the power 
spectrum and the squeezing spectrum have two peaks 
at the resonance frequencies u>k — v g(k 2 + k 2 ) 1 ^ 2 — 



1.545eV± [u 2 



t|/4] 1/2 as shown in Fig. 18b. The 



fiat regime in the Fig. 18b precisely and completely ex- 
plained the " anomaly " near k = of the excitation 
spectrum of exciton-polariton in a planar microcavity^I. 



APPENDIX B: THE TWO MODE SQUEEZED 
STATE FROM A MORE INTUITIVE VIEW 

In this appendix, we take a more intuitive way to show 
that no matter E(k) < 7^/2 or E(k) > 7s/2, the emitted 
photons are in a two mode squeezed state even off the res- 
onance. This conclusion is very robust and independent 
of any microscopic details. This appendix supplement 
the more formal discussions in section VI. 



7j/2, intersection point between the dashed blue line and the 
horizontal axis leads to the maximum squeezing which is a 
two mode squeezed state, but with non-zero squeezing angle 
shown as ± in Fig. 5. This horizontal axis leads to the the " 
anomaly " near k = of the excitation spectrum of exciton- 
polariton in a planar microcavity^. When E(k) > 7g/2, the 
intersection points between the dashed blue line and the solid 
gray lines lead to the maximum squeezing which is a two mode 
squeezed states with zero squeezing angle shown as ± in Fig. 5. 



The position and momentum ( quadrature phase ) op- 
erators of the output field can be more intuitively and 
straightforwardly defined as: 

X ± = Ag(a,)+Agt ( _ w)j 
Y± = (Bl) 
The Eq. ([HI]) and Eq. (55]) give the relation 

S k,X + ( UJ ) = = S ^( UJ ) and ^Y+M = 

St x (u>) = Sy(w) and the squeezing spectra: 



Sx{w) = 1 



S y (uj) = 1 + 



2nV d (kh 2 ^ 



2nV d {khl(e s + 2nV d {k)) 
Q 2 (uj)+j 2 E 2 (k) 



(B2) 



where = J 1 - [E 2 {k) - 7 |/4]. 

The two variance functions are defined as^: 

V±(u) =< X ± (u)Y ± {-ui)+Y ± (u)X ± (-w) > in /2 (B3) 

From Eqn l511 we can determine: 

v , ^ v f s 2 1% nV d {k)[u J 2 -E 2 {k) + ^] 
V+(td) = — V-\uj) = — - 



n 2 (u) +7 2 E 2 (fc) 



(B4) 
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From Eqn lB2[ we can find that: 



Sx{w)Sy(w) = 1 + 



^ 2 H(2nV^(fc) 7 g) 2 
[Q 2 (uj)+-/iE 2 {k)} 2 



> 1 



(B5) 



which seems suggest that the state may not be a min- 
imum uncertainty state. However, from Eqn lB5l and 
EanlBil we can see that AX±AY± = 1 + | V± | 2 even 
off the resonance which indicates it is a two mode 
squeezed state if we can choose quadrature phase prop- 
erly. From Eqn lB2l and Eqn lB4[ we can see that if we 
define: X± + iY± = (X ± + iY±)e ± ^ and choose squeez- 
ing parameter e = re t2 ^ such that: 



cosh 2r 



cos 20 



2(nVk(fc) 7 g) 2 

n 2 {Lo) + -fiE 2 {ky 



[n 2 (u;)+ 1 2 E 2 (k)} + (nV d (k)^) 2 



=(B6) 



Then we can show that Sx = AXl 



AX±AY± = 1 



-2r 



S Y = AY2 



(B7) 



which show that the emitted photons arc still in the two 
mode squeezed state in the basis of (X±,Y±) even off 
the resonance after making the rotation (f> in the original 
basis (X±,Y±) as shown in Fig. 6. It is easy to see Eqn lB6l 
is identical to Eqn[5Hl and Eqnl551 



APPENDIX C: THE RELATION BETWEEN THE 
CALCULATED QUANTITIES AND 
EXPERIMENT MEASURABLE QUANTITIES 

Note that the emitted photons at time t\ are described 
by the operator au{tx) in Eqn l431 so it is necessary to find 
the relation between afc(ii) and a2 u *(w) through the Eq. 
(|43p. The Fourier transformation of EqH3] leads to: 

v L_ 



k z 



where kl 



J 1 / 2 . So o2 ut ((jj) is a linear 



combination of a,^ kr Jti) and a,g _ s . r ^(^i)j both are at 
the output time t\. 
Then we can get: 

HuJ tl)a ^K)^ = <4° ut H a FV)}/^) 

= S 1 {k,cj)6(uj-uj')/D % (^(C2) 

where, in fact, only k' z = k z contributes and S\(k,uj) is 
given by Eqnl65l The total photon number at a given 



in-plane momentum k is : 

dw*I> s (u; fc )Si(£,w) =J2Si(k,u) (C3) 

k z 

which is nothing but the momentum distribution curve 
in Eqn l66l 

The total photon number at a given energy to = 
v g \Jk 2 + k 2 is : 

k k z k 

which is nothing but the energy distribution curve in 
EqnEl 

A putative single mode squeezing spectrum of ai ut (uj) 
can be achieved by coincident detection of two pho- 
tons both with in-plane momentum k, but one with 
z-direction momenta k z , the other with — k z shown in 
Fig. 16. As shown in sect V-A, in fact, the squeezed 
state is a two mode squeezed state between k and —k, 
so only the photons described by A^^oj) — [a|"*(o;) ± 

a™£(v)]/V2 = -^JDj:{uj + Ji)/2e iu)tl A° u * (v) are in the 
squeezed state where 

± (a ( _ Sifc:) (*i) +*(_*,-*;) ( C5 ) 
So the coincident detection experiment shown in Fig. 
15 between the photons with momentum (k, ±fcj) and 
(—k,±k z ) at the output time t\ can lead to the squeez- 
ing spectrum Sx(k, oj) in Eqn[5H Just as the Markov ap- 
proximation made in the section V, because the squeez- 
ing spectrum in Fig. 7 and Fig. 8 and power spectrum in 
Fig. 11 are narrowly peaked at to = and DAuj + fj,) 
is also a slowly varying function of to, so we can sim- 
ply set D j~(u) + n) ~ DAfj) as a constant prefactor at a 
given in-plane momentum k in the squeezing, power spec- 
trum, macroscopic super-radiance and correlation func- 
tion calculations in the section V. The photons described 
by A™* 1 (lo) is squeezed along the direction normal to the 

squeezed direction of the photons described by A~ ut (ui) 
in the phase space as shown in Fig. 5. 



APPENDIX D: GOLDEN RULE CALCULATION 
TO SECOND ORDER AND TWO MODE 
SQUEEZED PHOTONS 

This appendix was inspired by J.Ye's visit to Univ. of 
Texas at Austin. J. Ye acknowledge very inspiring dis- 
cussions with A. Macdonald, Qian Niu and Yao Wang. 
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There are also independent unpublished work by Jung- 
Jung Su and A. MacdonaldZ^. The input-output for- 
malism used in the main text is to solve a integral- 
differential equation under the standard Markov approx- 
imation which is assuming the coupling constant varies 
little over the characteristic energy which is the exciton 
energy fi ~ 1.54ei>, this approximation is valid here. So 
the coupling constant g k was treated non-perturbatively 
under the Markov approximation in the main text. To re- 
produce the results achieved in this formalism by Golden 
rule is very diffucult, although by pushing Golden rule 
to second order, one can see some features on two modes 
squeezing as shown in this appendix. 

1. First order Golden rule calculation 

1. The condensate at k — 

As shown in the main text, the initial state is taken to 
be \i) = \BEC)\0) with energy Ei = 0. Due to the in- 
plane momentum conservation, the final state contains 
one photon is |/) = \BEC)\lk z ) with energy Ef = uj kz — 
/i. The transition probability from the initial state to the 
final state can be calculated by the first order Golden 
rule: 

w * = yEi(«i 2 P/-^) 

= T N l0 (Dl) 

where H' is given by Eqn l25l and the exciton decay rate 
7o = \g 2 (k z )D below EqnEl 

From the Eqn IDll we can get the radiation rate from 
the condensate: 

P r ° d = N l0 n (D2) 

which is the same as Eqnl3"81 

It is easy to see that the first order Golden rule cal- 
culation is only a one photon process, but the photon 
coherent state discussed in section IV is a superposition 
of many photon number states, so it may be interesting 
to push the Golden rule calculation to second order. Note 
that the first order Golden calculation in££ just treated 
the exciton BEC as a two level atoms and did not use 
the BEC many body quasi-particle basis. 

2. The quasi-particles at k ^ 

As shown in the main text, the initial state is taken to 
be \i) = \BEC)\0) with energy E t = 0. Due to the in- 
plane momentum conservation, the final state contains 
one photon and one quasi-particle |/) = 7' -,\BEC)\lk) 

where k = (fc, k z ) with energy Ef = uj k — fi + u\k\. The 
transition probability from the initial state to the final 
state can be calculated by the first order Golden rule: 



n f 

= f E^WI^l4lo)l 2 l(^c|7_A-l^c)l 2 

x 8(u)k - M + u\k\) (D3) 

where H' is given by Eqnl2"5l 

Substituting the inverse of Eqn lI8l — U£7£ — -. 
into Eqn lD3l leads to: 

= — / dLo k Dj:{uj k )g 2 {k)vl8{uj k - n + u\k\\DA) 

where the photon density of states DAu)k) at a given in- 
plane momentum k is given in Eqn l44l Because u\k\ <C fi, 
then Eqn lD4l becomes: 

Wi = Y v bk (D5) 

where the exciton decay rate 7^ is given in Eqn l481 It is 

easy to see that the photon emission at k 7^ is indeed 
due to the quantum depletion in Eqn l23l 

From the Eqn lD4| we can get the radiation rate from 
the quasi-particles at a given in-plane momentum: con- 
densate: 

<(fc) = -jrtyj-H (D6) 

which cab be shown to be the same as Eqn[7D] to the 
linear order in 75. 

It is easy to see that the first order Golden rule calcu- 
lation is only a one photon process, but the photon two 
mode squeezing discussed in section V is a two photon 
entanglement between k and —k, so it may be interesting 
to push the Golden rule calculation to second order. Note 
that the first order Golden calculation in^i just treated 
the exciton BEC as a two level atoms and did not use 
the BEC many body quasi-particle basis. 

2. Second order Golden rule calculation 

The initial state is still taken to be \i) = \BEC)\0) 
with energy Ei — 0. Due to the in-plane momentum con- 
servation, the intermediate state |m) contains one pho- 
ton and one quasi-particle \m) — 7* J\BEC)\lk) with 

energy E m = u>k — ji + u\k\. AS shown in the last 
subsection, the matrix element between the initial \i) 
and the intermediate state |m) is V m i = (m\H'\i) = 
—ig{k)vx- There are two possible final states (1) two 
photons, but no quasi-particles: |/1) = \BEC)\l k )\l-k) 
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where k = (k,ki z ),—k — (— fe, &2z) with energy Ef 

1 1 ■ o;_fc — [i. The matrix element between the in- 
termediate and the final state is Vf m = {f^\H'\m) = 
ig{—k)uj: (2) two photons and two quasi-particles: |/2) = 



jl^\BEC)\l k )\l q ) where fc 



(k,k z ),q 



(q, <? z ) with 



energy Ef = u> k — fJ. + u\k\ + uj q — fi + u\q\. The matrix 
element between the intermediate and the final state is 
V fm = (f2\H'\m) = -ig{q)v^. 

The transition probability from the initial state to the 
final state can be calculated by the second order Golden 
rule calculation. The coefficient in front of the final state 
I/) ^ 



C 



(i) 
/ 







cf 



e„ 



E, 



}dt'(D7) 



When E m differes from E n and Ei, the contribution 
from second term leads to a rapid oscillation which does 
not lead to a transition probability growing with the time 
t, so the second term can be dropped. The tarnsition 
probability from the initial state i to the final state / is: 



where the sums is over all the possible final states, the 
energy conservation Ef = Ei, namely 2(uj kz — /x) = is 
enforced by the 5 function. 

From Eqns lD9l and IDlOl one can see that the emitted 
photons can only have a single momentum pinned at the 
chemical potential /i which is consistent with the results 
achieved in the section IV. In order to treat the pole 
structure in Eqn IDlOl one has to treat the condensate- 
photon by Heisenberg equation of motion as done in the 
section IV. 



b. The quasi-particles at k ^ 

1. The transition probability to the state |/1) = 
\BEC)\l k )\l. k ) 

The transition from the initial state \i) to the final state 
can be through two possible intermediate states 
\ml) = J_ % \BEC)\\ k ) and |m2) = 7 J|B£?C)|l_ fc ), so 
the transition probability to the final state after sum- 
ming over the two intermediate states is: 



Wi. 



2n 



\V fi 



^E 

S2Jm 



yfrri'TTi 



E, 



\ 2 5{E } - E^ (D8) 



W U1 



2tt 



1 



In the following, we will use this formula to calculate the 
transition to the state and |/2) respectively. 



The condensate at k = 



1. The transition probability to the state |/l)o = 
\BEC)\l kz )\l qi ),k z ^q z 

The transition from the initial state |i) to the final state 
|/l)o can be through two possible intermediate states 
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u\k\ 

2/t 



Lu k - M + u \k\ 



duj k D^{ui k )D^{2^ - uj k/ 
2u\k\ 



(uj k - - {u\k\) 2 



where the first and second sums are over all the possible 
ml} = \BEC)\l kx fand \m2) = \BEC)\l q J, so the tran- nnal states > the ener gy conservation E f = E l; namely, 



sition probability to the final state |/l)o after summing 
over the two intermediate states is: 



IF, 



(/l)o 



2tt 

T 



EE 



N 2 \g(k z )g(q z )f 



1 



1 



u kz - [i. 



8(tu kz - n + u>q z - u) = 



where the first and second sums are over all the possible 
final states with k z =/= q z , the energy conservation Ef — 

u — [i — is enforced by the S 



Ei, namely w kz 
function. 

2. The transition probability to the state |/2)o = 
\BEC)\2 kz ) 

The transition from the initial state \i) to the final 
state |/2)o can be through only one intermediate states 
|m) = \BEC)\l kz ), so the transition probability to the 
final state |/2)o is: 



^k — M + W- k — fi = 0, is enforced by the S function. . We 
can see that there are two resonances at to k — fi ± u\k\. 
We recovered the results achieved in Sec. VI-1 in the 
l^eak decay case E(k) > -f(k)/2 where the quasi-particles 
are well defined. This is expected because in the Golden 
(D9)rule calculations, we treat j(k) perturbatively, so can 
only discuss the weak decay case, the strong decay case 
is beyond the scope of the Golden rule calculation. Note 

that the prefactor upjj: — ^ d ffl is completely due to the 
exciton dipole-dipole interaction. 

In order to treat the pole structure near the two reso- 
nances, one has to treat the exciton-photon system self- 
consistently as done in the section V. 

2. The transition probability to the state |/2) = 



^(/^xE^v^ 



i 



Uk z - A* 



2 5{2{u kz 



-m)) 

(D10) 



71^^^)11^)11, 

The transition from the initial state \i) to the final state 
|/2) can be through two possible intermediate states 
|ml) = ^_j.\BEC)\l k ) and |m2) = jK\BEC)\l q ), so the 
transition probability to the final state |/1) after sum- 
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ming over the two intermediate states is: 

w (m = fEEw fc w«M J i — hr is 

11 kz q-_ UJ k -(l + u\k\ 

+ — ; 7-^\ 2 5(uJk - M + ll\k\ + U! q - fl + u\q\) 

= (D12) 

where the energy conservation Ef = Ei, namely, u>k — 
fx + u\k\ +Ljg — fx + u\q\ = is enforced by the S function, 
We conclude that there is no transition probability to the 
state |/2) in sharp contrast to the state 

3. Discussions 

Combining the calculations in the first and the second 
order Golden rule calculations in the last two subsections, 
we can write the complete Hilbert space of the exciton- 
photon system at a given inplane momentum k as: 

|*) = c Q \BEC)\0) +cs ll K\BEC)\l k ) 

+ c 2a |B£C)|l fe )|l-fc) 

+ c 2bl [^\BEC)\l k )\l q ) + --- (D13) 

where the • • • stand for multi-quasi-particles or multi- 
photon states. Eqn lD12l gives C2b = 0. It is the C2a 
term which leads to the two mode squeezing between 
k and —k. We can compare Eqn lD13l with the well 
known Wegner-Weiskoff theory of the spontaneous ra- 
diation of two level atoms where the Hilbert space is: 
I*) = c |e, 0) + C-.Ag, lg). Because the atom only has 
two levels: the excited state |e) and the ground state \g), 
so the atom+photon system can only has either no pho- 
ton or one photon state. In this case, the Golden rule 
calculation to first order can gives the correct decay rate. 
In the three level atom system with two photon cascades: 
|*) = c a \a,0)+c b k ~\b, lfc)+c c jr -|c, lg, the atom only 



has three levels: |a),|6),|c), so the atom+photon sys- 
tem can only has no photon, one photon and two photon 
states. In this case, the Golden rule calculation to sec- 
ond order is needed to give the correct decay rate. It 
is known there are no entanglements between the two 
photon states. 

While the exciton BEC system has infinite many body 
states with different number of quasi-particle excitations, 
so the exciton-photon system can has infinite many pho- 
ton states as shown in Eqn lD13l This is the crucial dif- 
ference between a quantum phase and a two levels or 
three levels system. So the Golden rule calculations are 
very limited in this case. To first order, Golden rule 
is just a one photon process, so it fails to capture the 
squeezed photons which are a two photon process. The 
two mode squeezed photons are the pairing of two pho- 
tons between k and — k just like pairing of two electrons 
between k and —k. In this appendix, we pushed the 
Golden rule to the second order which connects the initial 
and final state by two photons, we indeed find the signa- 
tures of the two photon squeezing. The formalism used 
in the main text is an input-output formalism which is to 
all order of the coupling constant, the only approxima- 
tion made is the Markov approximation which was also 
used in conventional Weisskopf-Wigner theory of sponta- 
neous emission^!. By looking at the squeezing spectrum 
Eqnl59l and the power spectrum Eqn l65l one can see the 
coupling constant 7g appear in both numerators and de- 
nominators which take the two photon squeezing into ac- 
count completely. The power spectrum which shows the 
macroscopic superradiance is also due to the two pho- 
ton pairing. Especially below EqnJSSl one can see the 
coupling constant 7^ is in the denominator in the MDC 

curve in the k — > limit which is in the macroscopic 
super-radiance case, so can not be achieved by the per- 
turbation Golden rule. One need to sum over infinite 
number of terms to reproduce the results achieved from 
the input-output formalism. 
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